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GULF  STREAM  GROUND  TRUTH  PROJECT: 
RESULTS  OF  THE  NRL  AIRBORNE  SENSORS 


I.  INTRODUCTION 

During  late  May  and  early  June  of  1976,  NASA  Wallops 
Flight  Center  (WFC)  conducted  a concerted  study  of  the  Gulf 
Stream.  The  experiment  combined  in  situ  measurement  of 
hydrographic  and  atmospheric  parameters  with  those  made  by 
a large  variety  of  airborne  remote  sensors.  The  sea  truth 
was  furnished  by  the  R.V.  Advance  II  (North  Carolina  State 
University).  NASA  Wallop's  C-54  aircraft  was  instrumented 
with  its  AAFE  SEASAT-A  altimeter,  NASA  Langley  Research 
Center's  AAFE  scatterometer  (RADSCAT) , Naval  Research  Labora- 
tory's (NRL's)  high  flight  radar,  wind-wave  radar,  a laser 
prof ilometer , and  a precision  radiation  thermometer  (PRT-5) . 


The  project  area  was  situated  due  east  of  WFC  in  the 
vicinity  of  the  Gulf  Stream.  This  location  was  selected 
for  the  expressed  purposes  of  (1)  testing  the  satellite 
altimeter's  ability  to  detect  geostrophic  current  systems, 

(2)  comparing  the  surface  temperature  signature  to  the 
actual  current  location,  and  (3)  checking  the  airborne 
altimeter's  sensitivity  to  rapid  changes  in  sea  state,  such 
spatial  inhomogeneities  in  sea  state  should  exist  there 
because  of  strong  wave-current  interactions  (Chen  and  Bey, 
1977).  Results  relevant  to  (1)  have  indicated  a positive 
capability  which  have  been  reported  by  Leitao,  et  al. 

(1977) . In  addition  to  these  goals,  the  experiment  afforded 
an  excellent  opportunity  to  resolve  some  fundamental 
questions  about  a prof ilometer ' s ability  to  obtain  correct 

Note:  Manuscript  submitted  April  26,  1978. 
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wave-frequency  spectra  as  well  as  information  on  the  wave 
field's  directional  properties.  NRL  was  also  responsible 
for  the  analysis  of  the  precision  radiation  temperature 
data.  The  results  presented  herein  are  only  those  derived 
by  the  NRL- sponsored  sensors  and  appropriate  comparisons 
are  made  with  those  of  other  devices. 

II.  EXPERIMENT  PLAN 

Because  the  aircraft  housed  several  instruments,  each 
with  its  own  particular  requirements,  a rather  complicated 
flight  plan  was  devised  by  the  participants.  Figure  II-l 
is  a depiction  of  the  general  scheme  and  the  sensors  to  be 
operated  at  each  altitude.  A 3048  m altitude  was  selected 
because  it  is  the  ceiling  for  flights  without  supplemental 
oxygen  and  it  is  necessary  to  fly  as  high  as  possible  for 
the  operation  of  altimeters  using  waveform  analysis  as  is 
shown  by  Walsh  (1977).  Star  patterns  of  different  geometries 
were  incorporated  at  the  1524  m and  305  m levels.  At 
1524  m,  the  RADSCAT  was  flown  upwind,  downwind  and  at  45° 
to  the  wind.  Figure  I I- 2 shows  the  proposed  pattern  at 
305  m where  the  NRL  wind-wave  radar  and  the  laser  profilom- 
eter  were  operated.  Two  sets  of  such  patterns  were  included, 
one  on  each  side  of  the  Gulf  Stream  front  with  continuous 
data  to  be  collected  by  all  designated  sensors  during  the 
3048  m transit. 

Three  missions  (5/27,  6/2,  6/4)  were  completed  during 
a two  week  period.  The  position  of  the  frontal  system  was 
supplied  prior  to  each  flight  by  the  Advance  II,  which 
continuously  maintained  stations  along  the  track  between 
the  coordinates  (36°28'N,  72°30"W)  and  (37°28.5'N,  72°30'W) 
during  this  time  interval.  Earlier  flights  included  an 
initial  systems  check-out  and  later  another  (5/26)  that  was 
4 scrubbed  due  to  a malfunction  in  the  Inertial  Navigation 

System  (INS) . 
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III.  INSTRUMENTATION 

NRL's  microwave  devices  mentioned  above  are  nadir- 
looking short  pulse  radars,  the  results  from  which  are 
derived  by  studying  the  return  waveforms  in  a manner 
appropriate  for  the  desired  product.  The  high-flight  radar 
provides  information  about  significant  wave  height  (SWH) 
from  the  slope  of  the  waveform's  ramp  as  described  by  Walsh 
(1974) . On  the  other  hand,  the  wind-wave  radar  infers 
surface  wind  speed  from  the  trailing  edge  of  the  waveform 
(Hammond,  et  al^.  , 1977).  Adjustments  to  the  results 
printed  in  that  paper  are  contained  in  Appendix  A. 

The  laser  profilometer  used  was  the  Spectra-Physics 
Geodolite  3A  similar  to  the  one  described  by  Ross,  et  al . 

(1970) . It  is  a continuous  output  (CW)  device  with  a 
selection  of  response  times  and  range  settings.  Because 
range  is  measured  by  comparing  the  phase  shift  between  the 
returning  amplitude-modulated  output  signal  and  a similar 
reference  signal,  range  ambiguities  exist  and  a phase  shift 
detection  subroutine  must  be  included  in  the  computer 
software  to  yield  a record  void  of  discontinuities.  These 
are  rather  easily  identified  because  their  magnitude  is 
similar  to  that  of  the  range  scale  being  used. 

Since  the  aircraft  (A/C)  is  moving  at  a high  velocity, 
the  effects  of  the  laser's  response  time  should  be  investi- 
gated so  as  to  understand  what  apparent  wavelengths  being 
profiled  will  be  diminished  in  amplitude.  Appendix  B 
includes  a derivation  of  the  cut-off  frequency  based  on  the 
definition  of  response  time  provided  by  the  operation 
manual.  Because  digitized  data  is  subject  to  Nyguist 
folding  and  the  associated  aliasing  of  energy,  an  additional 
calculation  using  the  appropriate  values  of  the  sampling 
j frequency  (f  ) and  cut-off  frequency  (f  ) is  included. 

| SC 

Additionally,  the  folding  frequency  (45  Hz)  is  high  enough 
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Fig.  Ill— 1 — Data  recording  system 
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Thus,  one  has  little  recourse  but  to  accept  this  limitation 
on  the  prof ilometer ' s accuracy. 

The  accelerometer  in  the  Wallops  INS  (a  Bendix  LTN  51) 
has  a specified  resolution  and  an  accuracy  of  0.001  g and 
0.01  g,  respectively.  A malfunction  in  this  INS  during  the 
5/26/76  flight  necessitated  the  borrowing  of  another  INS 
(LTN  72)  from  the  Patuxent  Naval  Air  Test  Center.  Unfortu- 
nately the  calibration  constant  for  that  system's  accele- 
rometer was  unknown.  However,  conversations  with  personnel 
from  Bendix  indicated  that  it  is  probably  the  same  instru- 
ment as  in  the  LTN  51  even  through  the  LTN  72  is  not 
normally  outfitted  with  one. 

Problems  with  the  Wallops  INS  continued  during  the 
6/2/76  mission  and  all  the  tracks  for  thfr  day  were  derived 
from  LORAN-A  fixes.  Large  disparities  between  the  INS  and 
LORAN-A  tracks  existed.  This  presents  a problem  since  a 
knowledge  of  0 (Figure  II-2)  is  essential  to  the  data  reduc- 
tion. Loran-derived  values  of  0 are  contained  in  Appendix 
C.  Data  on  6/4/76  offered  a higher  sea  state  and  was  quite 
adequate  for  the  wave  spectra  studies.  Also,  on  the  last 
flight,  excellent  agreement  existed  between  the  INS-derived 
wind  direction  and  both  the  ship's  surface  wind  and  dominant 
wave  directions. 

The  precision  radiation  therometer  is  a model  PRT-5 
manufactured  by  Barnes  Engineering  Company.  The  device 
operates  with  an  8-14  micron  window  (infrared) . The  field 
of  view  is  2°  which  at  3048  m produces  a 106  m spot  diameter. 
It  is  operative  in  the  temperature  domain  of  -20°C  to  75°C, 
but  has  a nonlinear  output  response  to  temperature.  Three 
range  intervals  are  available,  (-20°C,  15°C) , (10°C,  45°C) 

and  (40°C,  75°C)  with  two  full  scale  output  options,  0-1  V 
and  0-50  MV.  Our  application  fell  within  the  medium  range. 
Examination  of  Table  2 indicates  that  the  output  is 
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Table  2 — PRT— 5 Calibration 


TEMPERATURE 


10 

11 

12 

13 

14 


(°C) 


% FULL  SCALE  OUTPUT 


0.000 

2.430 

4.886 

7.366 

9.871 


15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


12.401 

14.956 

17.535 

20.140 

22.770 

25.425 

28.105 

30.810 

33.540 

36.295 

39.075 

41.881 

44.712 


INTERVAL 
USED  IN 
EQUATION  3 


28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


47.568 

50.449 
53.356 
56.288 
59.245 
62.228 
65.236 
68.269 
71.328 
74.412 
77.522 
80.656 
83.817 
87.003 
90.214 

93.450 
96.712 

100.000 
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particularly  nonlinear  near  the  high  temperature  end  of 
that  range.  The  data  in  Table  2 was  taken  from  the  operation 
manual.  For  this  reason  a temperature  subrange  (15°C, 

27°C) , was  selected  based  on  Figures  IV-1  and  IV-2  and  and 
a calibration  constant  was  calculated  to  be  0.371°C/%  full 
scale  over  that  interval.  The  0-1  V output  was  used  during 
the  Gulf  Stream  flights  and  since  the  recording  system 
registers  MV  values,  the  equation  for  temperature  (°C) 
becomes 

T = 15°  + (X  - 124.0)  (0.0371),  (3) 

where  X is  the  value  recorded. 

* Operation  in  cold  environments  (below  -20°C)  may  cause 

the  apparatus  to  provide  false  values.  A precision  thermis- 
tor Bolometer  mounted  in  a tightly  controlled  reference 
temperature  cavity  is  used  in  stabilizing  the  system.  If 
the  cavity  temperature  is  allowed  to  drop  below  the  design 
value,  incorrect  readings  will  result.  Also  the  wind  chill 
factor  can  elevate  the  minimum  ambient  operating  temperature. 

IV.  ENVIRONMENTAL  DATA 

When  attempting  to  draw  conclusions  pertaining  to  ocean 
wave  development,  environmental  history  is  of  primary  impor- 
tance. A large  number  of  factors  are  involved  in  the  physics 
of  this  phenomenon  and  the  technique  used  here  requires 
additional  information.  Consequently,  we  have  placed  much 
emphasis  on  the  collection  of  these  physical  parameters  and 
were  assisted  in  doing  so  by  several  different  agencies. 

The  data  gathered  included  surface  weather  charts  from  the 
National  Weather  Service  (NWS) , Gulf  Stream  frontal  analysis 
from  the  Naval  Oceanographic  Office  (NOO) , SWH  analysis  from 
Fleet  Numerical  Weather  Central  (FNWC)  and  additional 
supporting  data  taken  on  the  Advance  I I . Of  fundamental 
importance  are  the  surface  wind  speed,  direction,  duration 
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and  fetch,  air-sea  temperature  difference,  estimated  SWH  and 
the  direction  of  the  dominant  wave  component.  Other  factors 
such  as  the  curvature  of  the  wind  field  also  enter  into  SWH 
analyses.  Since  we  did  have  at  least  two  sources  of  SWH, 
the  ship  and  FNWC,  no  attempt  to  implement  a wave  height 
analyses  scheme  based  on  the  above  mentioned  data  was  made 
for  the  area  concerned.  Nonetheless,  that  additional  data 
is  included  in  the  figures  discussed  below. 

The  frontal  analysis  charts  indicate  that  the  flight 
patterns  were  flown  directly  over  the  Gulf  Stream  and  that 
the  location  of  the  current  was  stationary  in  the  test  area 
during  the  project  period.  The  experiment  area  is  indicated 
by  the  two  ship  survey  lines,  the  first  of  which  was  completed 
prior  to  the  flights.  The  in  situ  measurements  plotted 
along  the  ship's  track  are  in  excellent  agreement  with  the 
N00  data  derived  from  satellite  IR  determinations  (see 
Figures  IV-1  and  IV-2) . 

May  2 7 ; (Figures  IV- 3a  to  3d) 

The  general  weather  situation  at  the  beginning  of  the 
day  (0000  GMT)  was  as  follows:  the  experiment  area  was  under 
the  influence  of  a weak  high  pressure  cell  (1022  mb)  which 
was  the  eastern  portion  of  a high  pressure  ridge  extending 
from  the  Great  Lakes  region  southeastward  through  the  mid- 
Atlantic  states.  The  resulting  isobaric  pattern  was  in  a 
NE-SW  orientation  over  the  area  of  interest.  An  atmospheric 

r 

surface  wave  formed  at  that  time  approximately  1055  km  due 
east  of  the  area  along  a stationary  front  which  extended 
southwestward  through  central  Florida.  During  the  next  24 
hours,  this  wave  developed  into  a weak  low  pressure  system, 
intensified,  and  moved  rapidly  towards  the  northeast  away 
from  the  test  area.  As  this  low  migrated  into  the  North 
Atlantic,  the  high  pressure  ridge  extended  further  seaward 
from  the  Virginia  and  Carolina  coasts  causing  the  geostrophic 
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flow  over  the  test  area  to  rotate  to  a more  N to  NNE  direction 


by  2400  GMT.  With  the  weather  following  such  a pattern, 
fetch  and  duration  estimation  is  rather  subjective.  Indeed, 
swell  propagating  at  a large  angle  to  the  wind  direction 
(ship  report)  that  afternoon  definitely  confused  the  issue 
of  spectral  development  and  therefore  spectral  analysis  was 
not  routinely  made. 

June  2:  (Figures  IV-4a  to  4d) 

Initially,  a southbound  cold  front  extended  in  an  east- 
west  direction  just  north  of  the  40°  parallel.-  Also,  a 
surface  wave  which  had  formed  30  hours  earlier  and  had  been 
located  along  a stationary  front  (over  central  Missouri)  was 
now  centered  over  western  Virginia.  It  resided  along  the 
northern  extension  of  the  stationary  front  approximately  110 
km  south  of  the  advancing  cold  front.  At  this  time,  the  test 
area  was  experiencing  a southwesterly  flow  and  was  situated 
beneath  the  western  portion  of  a high  pressure  system  centered 
over  the  north  Atlantic  (31.5°N,  47°W) . About  1200  GMT  the 
cold  front  overtook  the  stationary  front  to  the  south  and 
became  stationary  itself.  It  now  extended  southwestward 
through  central  Maryland  and  Virginia  approximately  280  km 
north  of  the  experiment  area.  Thus,  a southwestward  flow 
continued  over  the  project  regions  until  the  cold  front  shifted 
southward  past  the  region  of  interest  sometime  between  0000 
GMT  and  0600  GMT  the  next  day.  The  result  of  these  atmo- 
spheric motions  lead  to  generally  steady  winds  with  a long 
fetch  oriented  towards  the  NE.  Therefore,  by  the  afternoon 
of  June  2,  the  seas  were  undoubtedly  full-developed  as 
definded  by  Pierson  and  Moskowitz  (1964).  Ship  reports  gave 
the  dominant  wave  heading  to  be  consistently  40  degrees. 

June  4 : (Figures  IV-5a  to  5d) 

After  the  above-mentioned  front  had  passed,  the  area 
remained  under  the  influence  of  a high  pressure  ridge 
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'ig.  IV— 4 — Meteorological  data  on  6/2/76  (Continued) 


r 


Fig.  IV— 6 — FNWC  significant 
wave  height  chart,  5/27/76 


Fig.  IV— 7 — FNWC  significant 
wave  height  chart,  6/2/76 


Fig.  IV— 8 — FNWC  significant 
wave  height  chart,  6/4/76 
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extending  in  a E-W  direction  over  the  Great  Lakes.  The 
associated  southwesterly  flow  continued  through  June  4.  The 
highest  SWH  during  the  project  occurred  on  this  day,  however 
the  low  wind  speeds  indicate  that  this  sea  state  is  due 
mostly  to  swell  generated  earlier  that  morning.  Rain  and 
clouds  did  hamper  the  laser's  operation  but  not  seriously 
and  the  wind-wave  radar's  waveguide  was  flooded  prohibiting 
its  operation.  The  ship  reported  the  dominant  wave  direction 
to  be  constant  at  all  stations  with  a heading  of  225°.  No 
weather  chart  for  1200  GMT  was  available. 


In  addition  to  hydrographic  and  atmospheric  data,  daily 
SWH  charts  are  available  and  are  presented  in  Figures  IV-6 
to  IV- 8.  The  contours  are  approximations  derived  from  ship 
and  buoy  reports  received  by  Fleet  Numerical  Weather  Central 
These  reports  are  subject  to  observer  errors  and  coding 
errors  during  transmission  and  are  carefully  edited  by  FWC. 
All  charts  refer  to  1200  GMT.  This  information  will  be 
included  in  the  comparisons  of  SWH  from  the  various  sensors 
and  the  ship  in  Section  VI. 


V. 


DATA  PROCESSING 


In  this  section  specifics  regarding  the  analysis  of 
profilometer  data  will  be  discussed.  Figure  V-l  is  a flow 
chart  of  the  algorithm.  The  analyses  of  other  sensors  used 
during  the  Gulf  Stream  Experiment  are  more  or  less  standard 
and  are  described  in  other  manuscripts.  The  power  spectra 
SL,  SA/C»  and  SL'  will  be  defined  in  the  later  part  of  this 
section. 

The  first  step,  that  of  unpacking  the  data,  concerns 
the  retrieval  of  the  laser  and  aircraft  vertical  accele- 
ration records  from  the  data  tape.  The  format  of  the  tape 
is  not  simple  since  a number  of  quantities  are  sequentially 
recorded  and  a masking  technique  is  necessary  for  picking 
out  particular  data  sets.  The  manner  in  which  the  data  was 
recorded  is  as  follows.  The  basic  data  block  is  called  a 
file . At  the  beginning  of  each  file  is  a header  which 
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DATA  ANALYSIS  FLOW  CHART  - AIRSPEC  5 


REMOVE  PHASE  SHIFTS 
FROM  LASER  RECORD 


UNPACK  DATA  FILE  FROM  ORIGINAL  TAPE- 
SCALE  LASER  AND  ACCELEROMETER  VALUES 


APPLY  HIGH-PASS  DIGITAL  FILTER 


BREAK  DATA  RECORDS  INTO  SEGMENTS  AND 
REMOVE  THE  MEAN  FROM  EACH  ONE 


CALCULATE  POWER  SPECTRA  S.  AND  SA/C 
FOR  EACH  SEGMENT 


CALCULATE  SWH  FROM  S 


DOPPLER  SHIFT  S, 


AVERAGE  SPECTRAL  ESTIMATES 


CONVOLUTE  SPECTRA  WITH  A HAMMING  FUNCTION 


TRANSFORM  ACCELERATION  SPECTRUM 
TO  DISPLACEMENT  SPECTRUM 


SET  0th  COMPONENTS  = 0 


APPLY  SL  > SA/C  CONDITION 


SUBTRACT  SA/C  FROM  S 


AVERAGE  SPECTRA  FROM  RELATED  DATA  FILES 


(*  see  Appendix  F) 

Fig.  V— 1 


! 
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contains  the  file  number,  date,  time,  altitude,  and  radar 
sampling  scope  settings.  Each  file  contains  32  records 
which  in  turn  contains  one  INS  latitude  and  longitude  fix. 

A record  also  contains  256  sweeps . A sweep  is  composed  of 
one  sample  from  each  of  seven  analog  input  channels  and  160 
samples  from  the  radar  sampling  scope  which  define  the 
recorded  radar  return  waveform.  The  sweep  rate  is  approxi- 
mately 90  Hz  which  implies  a sample  interval  of  1/90 
seconds.  The  analog  channels  (see  Figure  III-])  are  INS 
aircraft  pitch,  INS  aircraft  roll,  INS  aircraft  vertical 
acceleration,  laser  prof ilometer , time  of  day,  PRT-5  and  a 
spare.  The  eighth  input  is  the  radar.  The  pulse  rate  of 
the  radar  is  90  KHz,  but  the  system  is  designed  with  an 
integrating  circuit  which  averages  about  6 return  pulses 
before  the  signal  is  feed  into  the  analog-to-digital  converter 
which  is  triggered  at  15  KHz.  The  magic  number  of  167 
samples  per  sweep  is  derived  by  the  division,  15  KHz/90Hz. 

Therefore,  the  recorded  radar  pulse  is  comprised  of  a 
sample  from  each  of  160  consecutive  waveforms  from  the  A/D 
converter.  Each  file  contains  8192  = 2^  sweeps  and  an 

seven-track  tape  can  hold  seven  files.  After  unpacking, 
the  values  of  each  quantity  are  transformed  from  millivolts 
to  physical  units  using  the  appropriate  calibration  constants. 

After  the  laser  phase  shifts  have  been  removed,  both 
the  accelerometer  and  the  laser  records  are  high-pass 
filtered.  The  purpose  is  two  fold.  First,  any  spectral 
analysis  should  include  a "prewhitening"  of  the  data  to 
remove  segments  of  components  whose  periods  are  longer  than 
the  data  segment  and  also  to  decrease  the  sidelobes  of  the 
transform  that  are  introduced  by  analyzing  data  intervals 
of  finite  length.  The  second  reason  is  to  remove  aircraft 
motion  contamination  which  will  be  discussed  later. 

The  numerical  filter  applied  was  the  so-called  "Martin 
Filter"  (Martin,  1957).  It  is  a symmetric,  non-recursive 
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filter  which  features  a sine  termination  to  the  gain 
function  and  a correction  which  insures  unity  gain  at 
t = f_/f  = 0 (f  = cut-off  frequency,  f = sample  fre- 
quency) . The  sine  termination  is  introduced  to  avoid  large 
oscillations  in  the  gain  function  due  to  a sharp  cut-off 
(Gibbs  phenomenon) . Since  the  filter  is  a low-pass  filter, 
the  high-pass  filter  weights,  W„,  are  derived  from  the  low- 
pass  weights,  Wr  . 

L 

1 - WL  (0) 

WH(n)  = (4) 

WH  = -WL*n)  n “ 

There  are  2N+1  weights  in  the  filter.  As  with  any  filter, 
once  N exceeds  some  value,  it  becomes  more  efficient  to 
apply  the  filter  in  the  frequency  domain  rather  than  the 
time  domain.  The  other  two  parameters  required  to  generate 
the  weights  are  f and  the  "slope  of  weights",  h,  which 
originates  in  the  sine  termination  scheme. 

The  domain  of  x is  (0,  0.5).  With  the  value  of  f and 

S 

the  size  of  the  data  block  used,  the  value  of  x is  very 
small  and  this  can  cause  problems  with  the  filter  gain 
function.  Several  factors  determine  how  the  filter  is  to 
be  designed.  To  begin,  the  flight  patterns  for  the  laser 
included  five  legs  and  during  each  leg  three  files  of  data 
were  taken.  In  order  to  have  an  acceptable  confidence 
interval,  a compromise  had  to  be  reached  between  the  length 
and  the  number  of  data  segments  to  be  analyzed  and  averaged. 
Since  the  filter  removes  2N  points  from  each  file  and 
because  the  FFT  requires  2^  points,  the  selection  was  22.8 
second  data  blocks  = 1/4  file  or  2^  points  per  analysis. 
Therefore,  three  analyses  per  file  could  be  made,  and  nine 
spectral  estimates  per  leg  could  be  averaged.  The  22.8 
second  interval  was  chosen  because  it  was  felt  that  a 
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representative  sample  of  the  surface  would  be  obtained  over 
a distance  equal  to  22.8  S x 55  m/s  « 1250  m.  The  55  m/s 
figure  is  the  lowest  ground  speed  of  the  aircraft.  Thus,  a 
prewhitening  filter  would  require  f » 1/20  s ^ , (Davis, 

C j 

1974).  Since  f =90  Hz,  x = 5.6  x 10  . No  filter  with 

O 

N £ 1024  points  could  be  obtained  for  such  a x that  had  an 

actual  cut-off  anywhere  near  20  seconds.  The  alternatives 

are  to  increase  f , decrease  f , or  both. 

c s 

At  this  point,  aircraft  motion  plays  a part  in  filter 
design.  Originally,  the  accelerometer  was  to  be  used  as 
the  primary  method  for  aircraft  motion  removal.  However, 
most  aircraft  energy  has  periodicities  above  10  seconds  and 
would  be  concentrated  in  the  lower  harmonics  of  our  trans- 
form. These  are  the  least  accurately  estimated.  The  error 
in  these  harmonics  is  increased  when  the  power  spectrum  of 
the  acceleration  record  is  transformed  to  one  of  aircraft 
displacement.  This  transformation  is  accomplished  by 

4 

dividing  each  component  by  oj  , oj  is  the  circular  frequency. 
Therefore,  a higher  cut-off  frequency  than  0.05  s ^ was 
advisable.  It  should  also  be  remembered  that  the  accele- 
rometer was  not  of  the  highest  quality  in  the  first  place. 

Why  not  reduce  f by  digitizing  at  a lower  rate?  Many 
laser-prof ilometer  derived  ocean  wave  spectra  available  in 
the  literature  curiously  deviate  from  what  is  expected 
theoretically  and  from  other  empirical  methods  (Parsons  and 
Goodman,  1975) . Therefore,  we  opted  for  the  highest  sampling 
resolution  available  which  was  f = 90  Hz  rather  than  going 
to  a larger  spacing.  At  100  m/s  aircraft  speed,  a sample 
every  meter  would  be  obtained. 

With  these  considerations,  a trial  and  error  approach 
was  made  to  find  the  best  filter.  The  compromise  was 
between  the  10%  attenuation  point,  the  shape  of  the  filter 
ramp,  and  the  size  of  the  sidelobes.  Figure  V-2  shows  the 
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filter  gain  function.  The  10%  attenuation  point  is  at 
approximately  6 s.  The  filter  parameters  are  t = 0.0012346, 
h = 0.0001,  and  N = 512.  This  filter  requires  a great  deal 
of  timq  to  convolute  with  a data  file.  Actually,  the 
reduction  of  f by  a third  or  a fourth  would  sharpen  up  the 

S 

filter  response  greatly,  and  would  make  the  filter  more 

sensitive  to  the  value  of  f selected. 

c 

After  the  Fast  Fourier  Transform  (FFT)  and  application 
of  the  Hamming  convolution  to  each  spectrum,  the  above 
mentioned  transformation  is  carried  out  on  the  acceleration 
spectrum  to  produce  the  aircraft  displacement  spectra,  SA^C- 
Through  the  same  FFT  and  Hamming  function  the  power  spectra, 

S_ , derived  from  the  laser  profilometer  measurements  are 
also  evaluated.  It  should  be  noted  that  strong  objections 
may  be  raised  about  applying  a spectral  window  in  the  fre- 
quency domain  instead  of  the  time  domain  (see  Davis,  1974). 

In  this  instance,  the  filter  applied  to  the  data  was  probably 
adequate  for  reducing  the  sidelobes.  The  problem  is  further 
alleviated  by  the  fact  that  the  wave  spectra  are  not  char- 
acterized by  sharp  peaks  but  are  rather  smooth  with  gradual 
slopes.  Since  the  means  were  removed  from  each  data  segment 
prior  to  the  FFT,  no  energy  should  reside  in  the  zeroth 
harmonic.  However,  because  of  the  problems  with  the  filter 
ramp  and  the  relatively  large  spacing  between  adjacent 
harmonics  in  this  portion  of  the  spectrum,  the  zeroth  harmonic 
usually  did  have  a finite  value  after  the  Hamming  function 
convolution.  It  was  arbitrarily  set  to  zero. 


Another  condition  was  implemented  as  well.  No  component 
of  S A^c  should  exceed  SL*  Ideally  SL  should  always  be  equal 
to  or  greater  than  SA^C  because  it  is  the  sum  of  two  uncore- 
lated  motions.  If  SA^C  did  exceed  SL,  they  were  considered 
equal  and  due  to  the  aircraft  and  was  set  equal  to  zero. 
All  other  components  of  SA^C  were  subtracted  from  those  of 
SL  and  residual  spectrum  is  defined  as  SL' . The  significant 
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wave  height,  SWH,  was  calculated  from  the  result  given  in 
Neumann  and  Pierson  (1966) . 


H 


1/3 


- 


(5) 


where  cr  is  the  total  variance  or  power  in  the  residual 
spectrum,  S^. 

In  order  to  obtain  the  true  wave  spectrum,  must  be 
doppler  shifted  in  some  manner.  No  angular  spreading 
function  was  incorporated  here.  Instead,  two  approaches 
were  used.  The  first  approach  is  to  assume  that  all  waves 
travel  in  the  direction  of  the  ground  track,  and  V'  = +V 
depending  on  whether  the  track  was  actually  with  or  against 
the  waves.  The  second  was  to  assume  that  all  waves  traveled 
with  the  wind.  Then  the  aircraft  velocity  relative  to  the 
waves  would  be  V'  = V cos  0.  The  angle,  0,  is  defined  in 
Figure  II-2.  The  transformations  are  given  by 


= - 2_  + 
2V’  ‘ 


and 


\W 


gu,a 

V' 


4V'u 


1/2 


w = w 


(6) 


(7) 


where  "t"  and  "a"  subscripts  imply  "true"  and  "apparent" 
respectively. 

Finally,  the  values  of  all  nine  spectral  estimates  in 
each  leg  of  the  star  patterns  are  averaged  for  that  particu- 
lar leg.  To  determine  the  0.90  confidence  interval,  the 
result  found  in  Bendat  and  Piersol  (1971)  is  used.  Each 
spectral  estimate  affords  2 degrees  of  freedom  (Bath,  1974). 
Therefore  the  total  number  of  degrees  of  freedom  is  a 18. 


n 


„ a 
n / 2 


G(f) 

G(f) 


n 


(8) 


n,  (l-j) 


where  n = number  of  degrees  of  freedom,  a is  a parameter 


28 


determining  the  confidence  interval  desired,  (1-a) , G(f)  is 

A 2 

the  estimated  spectrum,  G(f)  is  the  true  spectrum,  and  x is 

the  chi-square  function.  Equation  (8)  yields  (0.62,  1.9) 
for  the  90%  interval. 

VI.  PRESENTATION  AND  DISCUSSION  OF  RESULTS 

A.  Wind  and 

Table  3 contains  all  of  the  results  for  wind  and 
significant  wave  height.  As  for  the  wind-wave  radar,  due  to 
the  malfunction  of  a transmitting  cable  on  5/27/76  and  clogged 
horn  by  rain  on  6/4/76  only  the  data  taken  on  6/2/76  was 
consistently  good  during  the  305  m patterns.  For  this 
reason,  only  one  data  point  was  obtained.  It  should  be 
noted  that  the  algorithms  for  the  wind-wave  radar  relates 
the  wind  speed  to  a 12.5  m reference  level,  while  the  ship's 
anemometer  is  positioned  24.5  m above  mean  sea  level.  The 
elevation  corresponding  to  the  NWS  values  is  not  known  to 
us.  The  aircraft  speeds  are  given  for  the  3048  m transects 
in  order  that  the  distance  between  measurements  can  be 
calculated.  The  actual  positions  and  orientations  of  the 
patterns  and  legs  are  shown  in  Figures  VI-1  to  VI-4b.  All 
of  the  3050  m tracks  proceed  in  a southerly  direction  and 
therefore  SWH  values  near  the  end  of  the  track  relate  more 
closely  to  Pattern  2.  Initial  values  from  the  high  flight 
radar  (3050  m)  were  obtained  in  locations  closer  to  Pattern 
5 although  a time  lag  of  3 to  4 hours  may  exist  between 
those  measurements.  The  laser  results  do  show  some  spatial 
inhomogeneities  between  legs,  ie.  Pattern  2 on  5/27/76.  The 
laser  data  is  most  consistent  on  6/2/76.  Unfortunately, 
there  were  no  satellite  overflights  coincident  to  these 
measurements. 

B.  Ocean  Wave  Spectra 

The  data  selected  for  spectral  analysis  was  the  data 
obtained  on  6/4/76.  Figures  VI-5a  to  VI-14b  are  spectra 
from  the  two  star  patterns  flown  that  day.  Each  figure 
contains  two  plots,  the  ocean  wave  spectrum  and  the  aircraft 
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Table  3 — Wind  and  SWH  Results 


DATE  NRL  LASER  SHIP 


PATTERN-LEG 

RADAR 

PROFILOMETER 

REPORT 

NWS/FWC 

COMMENTS 

m 

8 

8.7 

5-10 

Hl/3  (m) 

5/27/76 

13:21 

3048m 

2.45 

1.2-2. 4 

A/C  Speed* 

13:25 

3048m 

2.60 

88.3  m/s 

13:27 

3048m 

2.00 

<3  3048  m 

13:31 

3048m 

1.94 

15:06 

2-1 

1.23 

1.4 

2-2 

1.35 

15:24 

2-3 

1.80 

2-4 

1.26 

15:42 

2-5 

2.06 

6/2/76 

13:18 

3048m 

1.65 

1.2-2. 4 

A/C  Speed* 

13:23 

3048 

1.70 

85.5  m/s 

13:25 

3048m 

1.92 

'?  3048  m 

14:56 

2-1 

1.88 

1.8 

2-2 

1.72 

15:16 

2-3 

2.12 

1 

2-4 

1.78 

15:39 

2-5 

1.83 

3-UW 

1.90 

5-1 

1.98 

16:46 

5-2 

1.96 

5-3 

1.99 

17:05 

5-4 

1.90 

5-5 

1.89 

6/4/76 

13:18 

3048m 

3.27 

2. 4-3.0 

A/C  Speed* 

3048m 

3.26 

87.5  m/s 

3048m 

3.25 

@ 3048  m 

3048m 

3.29 

15:02 

2/1 

2.77 

2/2 

2.94 

15:28 

2/3 

2.84 

■ 

2/5 

2.64 

m 

15:51 

2/4 

2.36 

■ 

3-UW 

2.83 

l 

3-DW 

3.27 

16:48 

4-UW 

3.33 

K 

4-DW 

3.72 

Iv 

5-1 

3.34 

E 

17:14 

5-5 

3.32 

5-2 

2.67 

17:32 

5-3 

2.72 

5-4 

3.60 

1 

30 


Fig.  VI— 1 — Flight  tracks  at  3050  meters 


Fig.  VI— 2 — Star  patterns  on  5/27/76 


displacement  spectrum.  The  convention  for  "TAPE"  printed 
on  each  plot  is  data  tape  number — pattern  number  (Figure 
II-l) — leg  number.  The  plots  are  arranged  in  the  order  of 
increasing  aircraft  track  deviation  from  the  wave  direction. 
The  plots  are  log-log  covering  five  orders  of  magnitude  on 
the  ordinate  and  three  on  the  abscissa.  The  plots  of  the 
aircraft  motion  power  spectral  density  demonstrates  what 
was  discussed  in  Section  V regarding  the  problems  of 
estimating  that  parameter  using  short  data  lengths  (ft:  22.5  s)  . 
Since  the  energy  resides  in  a frequency  interval  where 
spectral  resolution  is  relatively  coarse  and  where  the 
transformation  from  acceleration  to  displacement  is  particu- 
larly sensitive,  the  proper  elimination  of  the  side  lobes 
introduced  by  the  discrete  finite  Fourier  transform  is 
difficult.  For  this  reason,  those  curves  do  not  break 
downward  on  the  low  frequency  end.  This  problem  might  have 
been  avoided  had  the  spectral  window  been  applied  in  the 
time  domain  before  the  mean  was  removed.  Also,  components 
with  amplitudes  too  low  are  not  plotted.  When  comparing 
the  standard  deviation  of  the  aircraft  motion  to  the  signi- 
ficant wave  height  (both  values  supplied  on  each  plot) , 
large  differences  in  the  first  is  often  seen  between  data 
segments  even  though  estimates  usually  show  little 

variation.  This  fact  indicates  that  the  algorithm  handles 
this  rather  perplexing  situation  reasonably  well.  Its 
performance  would  certainly  improve  if  the  data  set  could 
allow  for  longer  data  segments  and  a sharper  filter  gain 
cutoff  as  argued  in  Section  V.  Figure  VI-15  is  provided  in 
order  to  demonstrate  the  effect  of  erroneous  estimation  of 
aircraft  motion  on  the  SWH.  Appendix  D gives  the  derivation 
and  definitions  for  this  plot. 

For  each  leg,  two  spectra  are  shown.  The  "a"  desig- 
nation corresponds  to  a doppler  shift  which  assumes  that 
the  waves  are  traveling  in  a direction  parallel  to  the 


flight  track,  i.e.  0=0°,  180°.  The  "b"  relates  to  the 
assumption  that  the  waves  propagate  parallel  to  the  wind 
direction.  The  direction  of  the  dominant  wave  component  as 
given  by  the  ship  report  was  perfectly  downwind  on  that 
day.  However,  if  the  Pierson-Moskowitz  model  is  used  to 
relate  SWH  to  wind  speed, 

Hl/3  * °-0212  u19.5  ' (9) 

the  wind  speed  during  the  time  of  flight  is  insufficient  to 
generate  the  seas  that  were  encountered.  Indeed,  the  winds 
had  gradually  subsided  from  higher  values  (ship  report:  10- 
12  m/s  winds,  4-6  m seas)  early  that  morning.  Therefore, 
the  seas  were  not  in  a purely  generation  phase.  Nonethe- 
less, the  spectra  from  upwind  and  downwind  tracks  do 
conform  to  a 5:1  slope  as  expected. 

Included  are  two  additional  spectra  (Figures  VI-16, 

17)  taken  from  upwind  and  downwind  flights  near  Newfoundland 
during  February-March,  1976.  The  seas  were  much  higher  and 
offered  a good  test  of  the  algorithm’s  performance  because 
of  high  wind  speeds  experienced  by  the  aircraft  (up  to  40 
knots).  As  can  be  seen,  the  5:1  slope  is  maintained. 

Also,  when  these  are  compared  with  the  Gulf  Stream  spectra 
(lower  sea  state)  from  similar  tracks,  the  plots  align 
perfectly  showing  that  the  high  frequency  components  do 
saturate  and  are  maintained  as  increases.  Spectral 

grdwth  does  proceed  in  the  manner  predicted  by  Phillips  and 
aptly  depicted  in  Kitaigorodskii  (1973) . It  should  be 
noted  that  line-profiled  data  taken  from  a directionally 
dispersive  system  will  produce  power  spectra  that  are 
distorted  to  some  degree.  The  more  directional  the  wave 
field  is,  the  smaller  the  distortion  (see  Appendix  E) . The 
spectra  shown  here  tend  to  exhibit  the  least  distortion 
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Fig.  VI— 6a  . Fig.  VI-6b 


Fig.  VI 


Fig.  VI— 9b 


when  assumption  "a"  is  implemented.  If  the  seas  had  been 
dominanted  by  swell,  assumption  "b"  would  have  provided  the 
most  realistic  results. 

C.  Ocean  Surface  Temperature 

Figures  VI-18a  thru  19f  show  the  plots  of  surface 
temperature  in  °C  for  the  days  5/27  and  6/2/76.  Clouds 
obscured  the  surface  on  6/4  making  the  analysis  impossible. 
Clouds  do  appear  in  the  data  but  are  generally  easy  to 
recognize  since  they  produce  large  and  rapid  depressions  in 
the  temperature  record  which  cannot  be  attributed  to  any 
real  surface  features.  Some  small  depressions  exist  which 
may  be  due  to  light  clouds  or  may  be  cross-sections  of  cold 
water  bodies  from  the  continental  shelf  or  slope  entrained 
by  the  Gulf  Stream.  Each  individual  plot  represents  a data 
file,  seven  of  which  were  taken  each  day  along  the  3050  m 
track.  Only  the  first  six  files  are  presented.  Apparently 
on  both  days  the  mission  was  commenced  just  south  of  the 
Gulf  Stream  front  since  the  temperatures  are  generally 
constant  and  no  gradients  of  the  scales  shown  in  Figures  IV- 
1 and  2 were  measured. 

This  conclusion  is  supported  by  the  results  provided  by 
Curtin,  et  al.  (1978)  . Each  figure  represents  a track  only 
7-8  km  long.  Most  temperatures  related  to  the  ocean  surface 
lie  in  the  range  between  20°-24°C  while  ground  truth  measure- 
ments show  16°-17°C  water  west  of  the  front  and  a temperature 
maxima  near  25°C.  Our  results  are  not  corrected  for  atmo- 
spheric degradation  which  would  make  our  values  somewhat 
low. 
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VII.  SUMMARY 


The  results  of  this  experiment  show  that  the  active 
microwave  sensors  (high-flight  radar,  wind-wave  radar) 
provide  consistent  and  accurate  estimates  of  significant 
wave  height  and  surface  wind  speed,  respectively.  It  should 
be  emphasized  that  the  atmospheric  conditions  existing 
during  this  project  were  relatively  calm  with  wind  speed 
below  12  m/s  and  the  sea  states  not  higher  than  4 meters. 
Measurements  of  significant  wave  height  by  the  laser  profilom- 
eter,  in  general,  correlate  quite  well  with  those  from  the 
high-flight  radar  although  comparisons  between  simultaneous 
laser  and  radar  measurements  are  not  possible  since  they 
operate  at  different  altitudes.  The  laser  estimates  of  SWH 
from  individual  legs  of  the  star-patterns  flown  indicate 
that  considerable  spatial  variability  in  sea  state  can  exist 
at  least  for  low  sea  state  conditions. 

Satisfactory  removal  of  aircraft  motion  from  the  profilom- 
eter  data  was  a major  effort  in  this  study.  The  approach 
was  to  filter  the  laser  record  so  as  to  remove  the  low 
frequency  components  which  were  not  coincident  with  the 
ocean  surface  signal,  then  incorporate  the  aircraft's  vertical 
acceleration  information  for  the  removal  of  high  freqeuncy 
contamination  that  did  overlie  the  true  signal  obtained  in 
the  profilometer  record.  Examination  of  individual  spectral 
estimates  of  the  filtered  accelerometer  data  exhibited  large 
fluctuations  in  this  quantity,  yet  the  resultant  significant 
wave  height  estimates  remained  stable.  In  addition,  two 
results  of  profilometer  data  collected  from  another  ground 
truth  mission  (Newfoundland)  are  provided  to  further  establish 
confidence  in  the  algorithm's  performance  on  data  sets 
having  much  more  severe  aircraft  motion  contamination  and 
higher  sea  states.  The  results  are  very  satisfactory.  An 
alternative  approach  to  aircraft  motion  removal  is  to  take 
longer  data  segments  which  would  allow  an  adequate  estimation 


49 


of  the  entire  aircraft  motion  spectrum  and  subtract  this 
from  the  spectrum  provided  by  the  prof iiometer . This  approach 
would  eliminate  the  necessity  of  applying  a numerical  filter, 
but  would  require  much  longer  tracks. 

The  effects  of  angular  dispersion  on  prof ilometer- 
derived  ocean  wave  spectra  can  also  introduce  serious  problems 
since  the  flight  track  may  not  be  parallel  to  the  direction 
of  a wave  component's  propagation.  Experimental  results 
show  that  the  previous  practice  of  assuming  that  all  waves 
propagate  parallel  to  the  wind  can  cause  more  extreme  spectral 
distortion  than  simply  assuming  that  all  waves  travel  in  a 
direction  coincident  with  the  flight  track.  This  is  especially 
true  when  the  angle  between  the  track  and  the  wind  exceeds 
30  degrees.  However,  this  observation  would  not  be  true  if 
the  wave  field  were  highly  direction , then  the  previous 
assumption  would  be  more  appropriate.  To  clarify  this  issue, 
a mathematical  simulation  has  been  included  which  models  the 
spectrum  derived  from  a line- profile  sample  of  a directional 
wave  field.  Several  mathematical  forms  for  the  directional 
wave  spectra  were  used  and  track  angles  of  0,  15,  30,  45, 
and  60  degrees  to  the  wind  were  considered.  The  result 
indicates  that  spectral  aliasing  is  not  serious  for  wind- 
driven  seas  if  the  flight  track  is  within  15  degrees  of  the 
wind  direction.  Otherwise,  the  distortion  become  more 
prevalent  as  that  angle  increases.  This  distortion  may 
indicate  a possible  technique  for  determining  the  functional 
form  of  the  angular  spreading  function.  By  taking  airborne 
profilometer  data  at  various  angles  to  the  wind  simultaneously 
with  a stationary  wave  profiler,  various  angle  spreading 
functions  could  be  applied  to  the  airborne  profile-derived 
spectra  so  as  to  shift  it  back  upon  the  spectra  from  the 
stationary  sensor. 

Finally,  the  passive  airborne  temperature  sensor,  the 
PRT-5,  produced  very  realistic  values  of  the  sea  surface 


temperature  even  though  atmospheric  corrections  were  not 
employed.  The  temperature  range  from  PRT-5  lies  within  the 
domain  of  values  supplied  by  in-situ  measurements. 
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APPENDIX  A 


DETERMINATION  OF  WIND  FROM  THE  WIND-WAVE  RADAR 

The  results  published  by  Hammond,  et  al  (1977)  show 
how  wind  speed  can  be  obtained  from  the  return  waveform  of 
the  nanosecond  radar.  However,  the  particular  cases  derived 
in  that  paper  relate  to  the  aircraft  altitude  of  150  m.  It 
was  therefore  necessary  to  recalibrate  the  algorithm  for  an 
altitude  of  305  m.  Figures  A-l  and  A-2  represent  the 
appropriate  curves.  Figure  A-l  shows  the  normalized  wave- 
forms that  correspond  to  various  mean  wind  speeds  12.5 
meters  above  mean  sea  level.  The  wind  is  related  to  the 
slope  of  the  trailing  end  of  the  waveform  and  Figure  A-2 
plots  wind  speed  vs  slope  of  the  trailing  end.  The  device 
is  particularly  sensitive  at  low  wind  speeds  as  were 
experienced  during  the  Gulf  Stream  flights. 
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APPENDIX  B 


EFFECT  OF  LASER  RESPONSE  TIME 

The  Spectra-Physics  Geodolite  3A  has  a selection  of 
three  response  times,  5,  10,  and  20  ms.  The  response  time, 
x,  is  defined  as  the  time  interval  for  the  range  output 
signal  to  get  to  90  percent  completion  in  response  to  a 
step  change  in  range.  Thus  there  is  a smoothing  effect  on 
the  measured  profile,  the  severity  of  which  is  determined 
by  the  aircraft  speed  and  the  response  time.  During  the 
Gulf  Stream  experiment,  the  response  time,  x,  was  20  ms. 

If  the  aircraft  speed  were  75  m/s,  then  the  smoothing 
effect  covered  1.5  m. 

It  is  desirable  to  find  a cut-off  frequency,  f , so 

G 

as  to  determine  the  power  gain  at  the  Nyquist  frequency, 
f-  = f /2.  From  the  above  information,  the  time  constant 

r S 

of  the  laser  circuit  may  be  found* 

0.9  = 1 - e"aT  (1) 

Using  x = 0.02s,  a equals  115.  The  voltage  gain.  A,  is 
given  by 

A - [1  + {u/a)2]~1/2  (2) 

where  w is  the  radian  frequency  of  an  input  component.  The 
standard  definition  of  f corresponds  to  a voltage 
gain  = 0.707. 
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0.707  = [1  + (2irf  /a)2]~1/2  (3) 

C 

Immediately  one  sees  that  a = 2irfc  and  f = 18.3  Hz. 
Therefore, 

A = [1  + (f/fc)]-1/2  (4) 

The  power  gain,  G,  is  of  interest  because  distance  is 
directly  proportional  to  the  output  voltage  and  power 
spectra  derived  from  the  variance  of  the  distance  measure- 
ments is  the  final  product.  In  decibels,  the  power  gain  is 

Gdb  = 10  • log  (A2) 


If  fj  = 45  Hz,  then 

Gdfc  - -10  log  ( 1 + jfiy)  = -8.5  db 

which  is  the  case  for  the  analysis  used  in  this  study.  The 
laser  does  serve  as  a low-pass  filter  but  8.5  db  decrease 
is  not  generally  considered  sufficient  for  eliminating 
spectral  aliasing.  Fortunately,  for  our  application,  no 
significant  energy  is  found  at  wavelengths  w (75  m/s)  • 

(1/45) s < 2 m.  However,  if  the  data  is  digitized  at  slower 
rates,  then  that  wavelength  increases.  For  ff  = 22.5  Hz, 
Gdb  = -4  db  and  the  corresponding  wavelength  would  be  « 3.3 
m.  For  f^  = 11.25  Hz,  G^  = -1.4  db  and  the  wavelength 
increases  to  6.5  m.  The  maximum  amplitude  of  a 6.5  m wave 
is  approximately  0.5  m,  therefore  significant  error  could 
be  introduced  during  low  sea  states  if  the  data  is  digitized 
too  coarsely.  All  of  these  factors  should  be  taken  into 
j consideration  when  a profilometer  is  used. 
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APPENDIX  C 


DETERMINATION  OF  THETA 

As  discussed  in  Section  VI,  the  data  has  been  analyzed 
in  two  fashions,  the  analysis  using  assumption  "b"  requires 
the  angle  between  the  ground  track  and  some  reference 
direction  be  known  in  order  to  properly  doppler  shift  the 
wave  spectra.  This  angle,  0,  (Figure  II-2)  can  be  defined 
and  evaluated  in  several  ways  depending  on  which  vector  is 
used  as  the  reference.  These  choices  include  the  surface 
wind  given  by  the  ship  or  by  weather  maps,  wind  derived 
from  the  aircraft  INS  data  and  dominant  wave  direction 
provided  by  the  ship. 

Using  the  wind  direction  as  a reference  direction 
could  prove  troublesome  since  the  dominant  wave  may  be 
swell  not  correlated  to  the  wind.  This  can  cause  a large 
error  in  the  doppler  shift.  Since  the  wind  direction  can 
change  rapidly  in  space  and  time  due  to  the  presence  of 
fronts,  using  the  wind  to  derive  0 is  not  advisable  in  such 
situations.  The  weather  patterns  during  the  Gulf  Stream 
experiment  included  several  fronts  as  discussed  in  Section 
IV.  Furthermore,  use  of  INS-derived  wind  would  be  subject 
to  error  if  the  wind  vector  rotates  with  elevation.  Despite 
these  objections,  wind  vectors  from  the  INS  data  were 
derived  for  each  leg  of  each  star  pattern  on  the  6/2/76  and 
6/4/76  flights.  The  INS  provided  the  aircraft  heading, 
ground  track,  true  air  speed,  ground  speed  and  drift  angle. 
The  vector  analysis  was  not  made  on  5/27/76  data  because 
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the  aircraft  true  air  speed  was  not  recorded  by  the  INS 
operator.  The  best  way  of  utilizing  these  vectors  is  to 
average  the  individual  estimates  over  the  five  legs  to 
obtain  a mean  wind  for  the  star  pattern.  The  individual 
vectors  showed  a great  deal  of  scatter,  particularly  in 
direction.  The  data  on  6/2/76  produced  extremely  scattered 
wind  directions  where  a maximum  difference  of  138°  within 
one  star  pattern  existed.  The  INS  was  not  operating  properly 
that  day  and  Loran-A  data  was  collected  as  an  alternative. 

The  individual  vectors  from  the  star  patterns  of  6/4/76 
were  much  more  consistent  for  both  star  patterns.  The  mean 
wind  vectors  from  the  two  star  patterns  had  the  headings 
214°  and  225°  while  the  ship  reported  wind  and  dominant 
wave  direction  was  225°.  Therefore  on  6/4/76,  all  values 
of  0 were  derived  using  225°  as  the  reference  heading.  For 
the  reasons  stated  above,  the  dominant  wave  direction  was 
used  in  defining  0 even  when  the  direction  differed  from 
the  wind  direction  as  it  did  on  5/27  40-45°  misalignment 

reported  by  the  ship) . Also,  ground  tracks  derived  from 
the  Loran  A information  were  used  on  6/2.  Table  C-l  presents 
the  final  values  of  0. 
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TABLE  C-l 


VALUES  OF  THETA 


Pattern/ 

Ground 

Dominant 

True 

Design 

Day 

Leg 

Track 

Wave  Heading 

0 

0 

5/27/76 

2/1 

169.4 

225 

124.4 

180 

2/2 

283.7 

225 

238.7 

300 

2/3 

41.6 

225 

356.6 

60 

2/4 

134.1 

225 

89.1 

150 

2/5 

11.2 

225 

326.2 

30 

5/1 

161.5 

220 

116.5 

180 

5/2 

281.1 

220 

236.1 

300 

5/3 

50.6 

220 

10.6 

60 

5/4 

141.6 

220 

101.6 

150 

5/5 

20.5 

220 

340.5 

30 

6/2/76 

2/1 

0 

40 

140 

180 

2/2 

94 

40 

234 

300 

2/3 

240 

40 

20 

60 

2/4 

338 

40 

119 

150 

2/5 

192 

40 

336 

30 

5/1 

64 

40 

204 

180 

5/2 

210 

40 

313 

300 

5/3 

261 

40 

41 

60 

5/4 

13 

40 

153 

150 

5/5 

258 

40 

38 

30 

6/4/76 

2/1 

268.0 

225 

223.0 

180 

2/2 

39.3 

225 

354.3 

300 

2/3 

160.8 

225 

115.8 

60 

2/4 

242.7 

225 

197.7 

150 

2/5 

134.8 

225 

89.8 

30 

5/1 

232.0 

225 

184.0 

180 

5/2 

345.0 

225 

300.0 

300 

5/3 

128.2 

225 

83.2 

60 

5/4 

209.5 

225 

164.5 

150 

5/5 

93.2 

225 

48.2 

30 
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APPENDIX  D 


EFFECT  OF  ERRONEOUS  ESTIMATION  OF  AIRCRAFT  MOTION  ON  H.^ 


In  this  appendix,  the  derivation  of  n in  Figure  VI-21 

is  given.  The  subscripts  "m"  and  "T"  stand  for  "measured" 

2 

and  "true"  respectively.  The  symbol,  a , is  the  variance 
and  the  subscript  "L"  implies  "laser  and  "A/C"  relates  to 
the  aircraft  motion. 


(H1/3>T  ’ 4 [°l  - 
(Hl/3>m  " 4 [°L  - 

<Hl/3>m  _ ft 

(h1^3)t  |^a2 


«4/c>t] 

<<=>.]  1/2 

- <<’j/c)m] 

^ °A/C  ^ T J 


1/2 


If  " * 


and  e = (0l/C)m/(0i/C)T'  then 


_<e  (n-l)  \ 
len-l  / 


1/2 


Thus  the  ratio  of  measured  to  true  SWH  is  a function  of  the 
relative  size  of  the  measured  aircraft  motion  to  the  true 
aircraft  motion  and  the  ratio  of  the  laser  measurement  to 
the  measured  aircraft  motion.  Figure  VI-21  presents  the 
family  of  curves  for  n = constant. 


APPENDIX  E 


EFFECTS  OF  WAVE  DIRECTIONALITY 
ON  LINE-PROFILE  DERIVED  SPECTRA 


We  will  simulate  analytically  the  frequency  spectrum 
that  would  be  obtained  by  using  surface  elevation  data  taken 
by  a profilometer  moving  in  a straight  line.  Angular 
dispersion  in  the  wavefield  will  affect  our  results.  The 
wave  number  of  any  wave  component  will  be  alaised  in  the 
manner  indicated  in  Figure  E-l.  Please  note  that  the  0 
defined  in  Figure  E-l  is  not  related  to  the  0 defined  in 
Figure  II-2  of  the  main  text  and  is  used  here  in  order  to  be 
consistent  with  most  of  the  literature  describing  angular 
dispersion  of  ocean  waves.  To  derive  a unique  quantitative 
result  is  not  possible  since  no  single  form  for  the  direc- 
tional wave  spectrum  has  been  proven  to  be  superior  to  other 
proposed  expressions  found  in  the  literature.  Therefore,  we 
will  select  a model  spectrum  from  the  literature  that  yields 
qualitative  answers  to  this  question. 

The  directional  wave  spectrum  is  generally  written  as 
the  product  of  a frequency-dependent  angular  spreading 
function,  H(a,0),  and  a frequency  dependent  function. 

$(o,0)  = H(o,0)<t>'  (a)  (1) 

This  form  has  been  suggested  by  many  investigators 
including  Lonquet-Higgins  (1962)  and  more  recently  Ou,  et 
al.  (1974),  but  insufficient  data  exists  to  really  specify 
H (o , 0 ) . More  often  a simplified  version,  h(0),  is  used.  On 


E-l 


the  assumption  that  wave  energy  should  be  distributed 

symmetrically  about  the  wind  vector,  h(0)  is  usually  chosen 

to  be  a cosine  function  raised  to  some  power  or  a summation 

of  such  terms.  Indeed,  Krylov,  et  al^  (1966)  has  indicated 
2 

h(8)  cos  0.  Therefore,  we  will  use  the  following  functions: 

-1 


h(0)  = B (n)  cos2n  0 = 


TT 


Isotropic 


(2/tt)  cos  0 
(8/3tt)  cos^  0 


n=0 

n=l 

n=2 


(2) 


where 


B (n)  = 


tt/2 


i-1 


cos 2n  0 d0 


■ltt/2 


B(n)  represents  a normalization  factor  which  ensures  the 
conservation  of  energy. 

The  equations  developed  below  are  general  in  that  0 ' 
(see  Figure  E-l)  may  assume  any  value  within  the  domain 
[ — tt/2 , tt/2].  The  relationship  between  <p'(a)  and  the  fre- 
quency spectrum,  <J>(a),  is  found  by  integrating  4(0,0) 
with  respect  to  0 and  <t>'(a)  = 4>(a)/a.  A similar  relation- 
ship holds  in  wave-number  space.  S' (k)  = S(k)/k. 


The  Pierson-Moskowitz  (P-M)  spectrum  will  be  used  for 


<t>  (o)  • 


M 


<Ha)  = 


20- 


exp 


(3) 


where  g is  the  gravitational  acceleration,  U is  the  surface 

_3 

wind  speed  at  19.5  meters  and  a and  B are  4.05  x 10  and 
0.74  respectively. 

Since  the  aliasing  problem  deals  with  spatial  rather 
than  temporal  frequencies,  4»(o)  is  transformed  to  S(k)  via 


E-2 


MMI 


r 


i 


the  dispersion  relationship,  a - gk. 


S(k) 


(4) 


The  directional  wave  number  spectrum  may  be  written  in  a 
form  analogous  to  (1)  . 


i,(k,e) 


(5) 


i 


From  Figure  E-l,  the  apparent  wave  number,  ka,  is  related  to 
the  true  wave  number  by 

k = k cos  y , (6) 

a 


and  therefore, 


S, (k  ) = S (k)  sec  y. 

a a 


(7) 


Substituting  (2),  (6),  and  (7)  into  (5), 


f (k  ,0)  = h(0)  • S'  (k  ) 

d a d d 


cos2n  (0'+y)  cos3y  exp 


— vcosy 


(8) 


where  A = oB(n)/2  and  v = 0g2/U4.  It  should  be  noted  that 
Q • + y = 0 for  all  cases  only  when  y is  defined  as  the  angle 
measured  from  the  flight  track  to  the  wave  ray.  To  obtain 
the  desired  result,  S, (k  ) , (8)  must  be  integrated  over  the 

coordinate  plane  from  0 = -it/2  to  tt/2.  Because  0'  = constant 
for  a particular  track,  d0  = dy.  Thus, 


E-3 


(9) 


Sa  (*a. 
a a 


+tt/2 

'0'  > - f W 9'' 

~tt/2 


Y)  ka  dy 


and  the  apparent  wave  number  spectrum  retains  an  additional 

dependence  on  0 ’ . At  this  point,  specific  values  of  n must 

be  substituted  into  (8)  in  order  to  determine  (9) . The 

results  involve  modified  Bessel  functions  1^  (C)  where 

C = 8(g/U2k  ) 2/2  and  the  integrals  are  given  by  Olver  (1964) 
a 

after  having  been  simplified  to  the  appropriate  forms.  It 
should  be  remembered  that  these  functions  are  even  when  j is 
an  even  integer  and  odd  when  j is  odd.  The  expressions  for 
the  apparent  wave-number  spectra  are 


ae 


-C 


4k 


3“  (Io  ~ h (?)) 


n=0 


S (k  ,0’)  = 
a a 


ae 


8k 


a 

-C 

3 . (2  + cos  2 0 1 ) • Iq  ( £ ) - 2(l+cos30')  •I1(0 


[' 


+ cos  20 ' • I2 (? ) 


ae 


-C 


24k 


D 


n=] 


(10) 


3 j 2 ( 3 + 2 cos  29')  IQU)  - (8  cos  29’ 


+ cos  40'  + 6) *1^  (C)  + 2 (cos  49'  + 2 cos  29') 


• I2(£)  - cos  40'*I3(c)J  n=2 , 


In  order  to  present  the  results  in  a form  that  generally 

appears  in  the  literature,  we  have  transformed  S (k  ,0') 

a a 

back  to  the  frequency  domain  and  nondimensionalized  those 
expressions  in  the  same  manner  as  Pierson  and  Moskowitz 
(1964)  . The  transformation  is  given  by 


E-4 


(11) 


*a<fa' 


0') 


ai  = 

u5 


4tt 


U' 


ftf 


Sa(ka< 


0') 


Figures  E-2  thru  E-6  show  the  apparent  frequency  spectra 

for  various  track  angles  relative  to  the  dominate  wave 

direction.  The  P-M  spectra  is  plotted  as  a reference  and 

the  isotropic  case  is  included  because  it  is  the  limiting 

case  as  directionality  decreases.  Several  trends  can  be 

seen:  the  spectral  peak  shifts  to  lower  frequencies  for 

larger  values  of  0 • ; the  peak  broadens  and  decreases  mono- 

tonically  as  0*  increases;  the  slope  of  the  high  frequency 

portion  of  the  apparent  spectrum  does  not  dramatically 

deviate  from  that  of  the  P-M  spectrum  until  large  values  of 

2 4 

0'  are  reached;  the  difference  between  the  cos  0 and  cos  0 
directionality  models  is  not  particularly  strong  until  0' 
exceeds  45  degrees.  The  spectral  distortion  is  not  severe 
for  track  angles  within  15  degrees  of  the  wind  vector. 

These  indicate  that  airborne  profilometer  data  is  still 
useful  for  wind-wave  generation  studies  as  long  as  close 
attention  is  given  to  the  track  angle  relative  to  the 
dominate  surface  wave  direction.  Also,  the  proposition  of 
using  spectra  from  various  track  angles  to  infer  h(0)  does 
not  appear  promising  because  the  spectra  do  not  show  a 
substantial  enough  variation  with  h(0). 
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AIRSPEC  5 


Ha  00  PRlNl , /ll O.MCC.A IN,S 
Ha  00  tQUIP,) «*»,MJ . MP 
Mi  U0  E3JlP.b»**.MY, AO 
Mb  po  e3 J 1 F . 1 0 = • • , h I , PR 
MS  on  r 1 \ , X • I 

PROi.RAM  AHSPfcC? 

I'  I Mr  MS  ION  YlM.92).Z(8L9>'>,ASQi'NU<409/.2> 

UlM.NSION  LQKESFIS)  #FHFACTOR(b> 

DIM:  NS  I ON  PL T ARRAY  11028), I M < 2 ) , 0 < t>l94  > 

dimension  ►nuht<409/,2> 

dim  NSION  f RPQSPPC140V7.2) , HbARF ( 2 > 

Dll.  NS  I ON  91102S) 

Di Mr  NS  JON  LPHKIO), LSI  12),  INDl?>.l)FNC(30l.2> 

0 I Mi-  NSION  r 1 1 T<2S/).WHGT<2!>7),  1FILT13) 

U i Mr  NS  I ON  UPON<301).LABOVF<2),LBtLOM(2) 

DIM:  NS  I ON  NGSPECUS) , EMAX < 2 ) , FBAR1 2 ) , SpbCAVGNl 2 > 

COMhON/A/Y, Z 

C(  MM3M/ROrAPPrt/LFILt<  /)  .KMUF.  I ► ILE.MFILE 
n(iMMDN/llN>'.K«/RAD(  ion  ) , CHAN<7) , AINSS18) . iNCOCt  (8) , NFlLfc, 

INLAY  ,MonT  m.  NYR.NHR,  MIN.  ALT,  UOuT.NHSCOPE.  NySCCPP.  VSCAl-E# I SWEEP 
COM'-ON/Sf  /LPrN,D£LT#PLTSPACE 
COM:  3 N/3/  y nr 

COM: ON/B/ J( PNOn, 2 ),HniHlF< 4097.2) 

ECU  VAl.EN>  (O.FnUMT),  (FILT.WElGT),  (U.)  NUHTF) 

TYP  REAL  KURT0SIS,LA8UVE,LfifcL UM 

TYPe  INTE3PP  PRINT, PHAS.Pl.TON, THETA 

DaTa  (IFI.T  s Hh  LOW.BH  HIGH.BH  FU) 

Da  T a ( (FpF  AC  TOR  f HP )» HP  b 1,?)i  1 . , U . 5. 0 • 33, 0 , 1 S, 0 , 0 5) 

C 

C NOS  » NOMjtR  OK  SAMPLES 

C IMj)  i NJMHER  OF  POINTS  IN  FAST  FOURIER  TRANSFORM 

C IP(F)  a PJwfcR  OF  2 FOR  1 M ( 1 ) 

C ISfcl  = NU4HER  OF  DAI  A SAMPLES  FROM  WHICH  ONE  DATA  »UlNT  IS 

C St .FC  TEn  FOR  SPECTRAL  ANALYSIS 

C XMAXVL  B MAXI  MUM  VALUE  UF  FHFUUENCV  FOR  PLOTTING  S»ECTKa 

C L’EL T = S»MPLP  INTERVAL  CALCULATED  FROM  THE  SFtED  0-'  THE  TAPE 

C FECOPl'ili 

C V B AIRCRAFT  GROUND  SPEED  IN  KNUTS,  IF  = 0»  NO  D0PsLER  SHIFT,  IT 

C SH0UL3  BE  NOTED  THAI  CONST  CONTAINS  A DIVISION  BT  V. 

C MB  MIND  SPEED  IN  METERS/SEC 

C IFIlDATbI,  FILTERED  Data  IS  RECURUED  on  OuTPIT  UK: 

C COT  b CUT*PFF  FREQUENOY/SAHPLlNU  FREQUENCY 

C H b F | LTt  T SLOPE  OF  MEIGHTS 

C NUMi-T  B NJMRER  OF  FILTER  HEIGHTS,  2N*1bNUPBEH  qT  »TS  In  FILTER, 

C Lt'Mi-1  a L jW  PASS  FILTER  « 1,  HIGH  PASS  FILTER  ■ Z. 

C NO  PASS  FILTER  « 3 

C Na Vi>L  a Same  AS  NOS 

C LPHi  b TITLE  OF  THE  RUN 

C FPELSPECbaVG  SPECTRUM  OF  TOTAL  NUMBER  OF  SPEORA  COMPUTED 

C PLlON«  C.  SPECTRA  ARF  CALCULATED.  LISTED  IN  OUTPUT  AS  PLOT. 

C SMHmNLVbI ,i , SWH  IS  THE  ONLY  UUANTITY  DERIVED  FROM  FFT  ESTIMATES 

C PHDiSTbO.],  PLOT  HAVE  HEIGHT  distribution 

C NC'MtVPLl*!.  NO  HISTUGRAhS.  LISTED  IN  OUTPUT  AS  HAVE  PL, 

C NUMhISTb  MLASPR  HISTOGRAM  PLUTTFU),  b2 ( ACCEL ERAT I 3N  PLOTTED  ALSO) 

C NIB] , plot  f reospec 

C N2  S 1 IN;nRPORATES  ANGULAR  SPREADING  INTO  DCPPlEK  SHlFl 

C N.Abi,  a/C  ilSPL  SPECTRUM  IS  REMUVEU  FROM  LASER  SPEDTRUM 

C IpAoSE  B } B NO  MACHINE  PAUSP.  ■ ANY  OTHER  IMEGEH  a pAuSE,  660U 

t PRInT  b 1j  PRInT  b V.  NOT  Tu  PRINT  b l 


F-l 


C FHAS  * SOT  TO  REMOVE  PHASE  SHIFT  ■ 0.  TO  REMOVE  PHASE  SHIFT  • 1 

C 1 F L C;  T S P * ).  SPECTRA  ( ASU2NU ) ARE  PLOTTED 

C NKEi  SAP  * NUMBER  OF  FILLS  TO  SKIP  ON  OUTPUT  TAPE  FILES 

C Mp£\  « COJED  NUM0ER  FOR  COLOR  PENS  OF  THE  PLCTTER 

C NkEt  a SUHBER  OF  RELOHDS  PEH  FILE  ON  INPUT  TAPE 

C LI  = INPUT  TAPE 

C LO  = OUTPJT  TAPE 

C PlTSPACE*jARAMETER  IN  PL0T2,  jT  EQUALS  INCHES  PER  3ELT . 

C MlRhW  ■ Td  REWIND  INPUT  TAPE  ■ 1.  NOT  TO  REWIND  IN»UT  T*PE  * 0 

C THETA  ■ ANGLE  IN  DEUREE  BETWEEN  GROUND  track  and  Mi aN  WiNU 

C DIRECT  ION 

C OP  = LASER  range  SCALE 

C Df LASER  s CONVERSION  COEFFICIENT  FOH  LASER  CCHRESPONDjNu  TO  OP  . 

C IT  IS  SET  NEGATIVE  TO  INVERT  THE  LASER  PROFILE, 

C DF  = CONVERSION  COEFFICIENT  FOR  ACCELEROMETER 
C NSPrCAVG  • NUMBER  OF  SPECTRA  TO  BE  AVERAGED  TOGETHER  PER  FILE 

C NGSFEC  ■ NUMBER  OF  SPEC1RA  IN  SEQUENTIAL  NUMEfcR  WHICH  AWE  NOT  TO 

C -E  AVERAGED 

C TSPFCAVG=TOTAL  NUMBER  OF  LASER  SPECTRA  DERIVED  TO  VjELO  FREOSPEC. 

C VSHIP  > SHIP  SPEED  IN  FEET/SEC 

C NF  = NUMBER  OF  FILES  10  BE  UNPACKED 

C KlC  s NUMTFr  of  THF  laser  CHANNEL 
C NaC  a NUMBER  OF  THE  ACCELERATION  CHANNEL 

C LFliE  a PARAMETER  DETERMINING  HHJCH  FILES  ARE  UNPACKED 

C KFltE  s NJMHER  OF  FILE  INPUT  TAPE  IS  PUSlTlCNEU  T3  READ 

C ISWrEP  a NUMBER  OF  SWEEPS  PER  RkCURD 

c NFLfG  a parameter  determining  no.  uata  pts. /sweep  to  be  unpacked 

C NsE7S=T0T4-  NUMBER  UF  DATA  sets  TO  BF  ANALYZED  SEPARATELY. 

C RES  I ART  = 1,00  INDICATES  THAT  COMPUTATIONS  OK  A DATA  SE i 

C ( WHICH  hay  be  COMPOSFO  OF  SEVERAL  F I LES ) ARE  F IN  j SHED  AND 

C T NAT  REINITIALIZATION  OF  PARAMETERS  FOR  NEXT  DA i A SET 

C MJST  BE  DONE. 

C 

9001  format  (SjR.rio.z.FS.o.jiFifi.s.Ais.Fio.S) 

90U2  FuRFAT  (1JA8) 
vous  format  <iii5> 

VOuA  FORMAT  ( // » ?u X » *THE  CONTROL  CAROS  FOR  THE  Nt  XT  F l _E  ARE ,.*,//, lnX 

• .•On  CONTROL  CARO  l -•,//, 5X. •NOS«JX .* IM < 1 )*JX. *1 H( 2 ) *6X. *DELT*6X . 
••V*tX,»CUr«8X,«H*AX»«NUMpT*3X,*LOWHI«6X.*NF»3X,«NANGL*AX*«DF*»/»SI 
«B,F' 0.7.F S.2.2F10.5#«I8,n.0,5,//,lOX.*ON  CONTROL  CARD  2 . LABEL  TN 

• 80  COLUMRS»//2X.10A8//10X*ON  CONTROL  CARD  .1  -*//6X*Nl*oX*N2*3X*<»A 
•USE»3X»pRINT«3X.pHASE«4X»pL0T  NHECSKP  WaVE  pL»5X*PEN*2X.Sp  AVG*,1X 
*.»lt L0TSPa//1118//lUX*0N  CONTROL  CARD  4 -*//3X*NUM  r£C*JX, 

• 7HLU  HPT,4X.*LU  0U1*;»X,12HPL0T2  PT  I NT  , t X , 4HMTRW , Bx , ZHUP . 5X  . 

• •1H(-TA*.  1X.»DF  LASER*//,  I1U,  2 1 1 0 , F 15 . 8 . 1 1 D , FlO  . 0 , 1 1 0 , Fl  4 . 

•3.//) 

90u6  FORMAT  <»J  DATA  FILE  CONTAINED  *I5»«  SAMFLES*/) 

90 U / Format  (•  *,5X,*MAX  VALUE  a *,FB.2,5X,*MIn  VALUE  s **18.2) 

90ub  format  n ri .rrx.iahpower  spectrum//zox.i2hthis  run  is  ixoas./.isx, 

• 16Hm UMBER  TIF  OBS  ■# In,50X,14HDATA  FACTOR  a . F 1 0 . 3 , / • 15X , 16HNUHBER 

• (|F  -AGS  al6.50X»»CUT»9X,«a  *F9,4»/l5X,*DELTA  T l ME **»X , •« *F14 , /,  4PX 

• •SLOPE  »7X,  *s  *F9.4./1!>X«A/C  VELUC  1 T Y*3X , *B«F1 1 . 4 . 4 5X*WE  i GHTS*5X , ♦ • 

• •I-.,/15XtVARlANCEU)  ■*Fll.4,45X*H  1/3  (1)  a *F  V , 4 , /1SX«  VA  9 I 

• AFCt-(Z)  a.FU.4,45X«STD  DEV  A/C  « «F9 . 4 , /4  VX , AH , * PASS  FILTERED 

• */  ) 

9!)U9  Format  <•  APPARENT  variance  apparent  CIRC  TRUE  C 1 PC 

• TRUE  TRUE  FREQ  TRUE  TRUE  WAVE*/*  FRfc Q*1 9X , «F Rfc 

•Q  SmEC  FREQ  FREQ  SPFC  FREQ  SPECTRUM  WAVE 

• Ni.RBER«/A3X,»hZ*.  BOX,  •numb  spectrum*//! 

9010  format(ihi,4x,*freqimz)  vahiance  fheq  sfec  :i«c  freq  circ 


F-2 


n o o c.  o 


T 


• fri-o  spfc :•.//» 

Mill  FORmAT(4E12.4,E1T.«) 

90l2  FORMAT  (14".//,5X,«IST(«. II,*)  ■ • . 1 5, 5X , • I NC ( • , 1 1 . • > ■ *»I5,/.5V. 
••OEwTROlU  * »Fl o . 4 , /5X.STD  NEV  ■ *F10 . 4 . /SX«SKEPNisS  ■ *F1U,4,/5X 
••KUPTOSlS  * *F  1 0 . 4 ) 

014  FORMAT  (3IS.  F 1U.B. 15.F5.U. IS.Fb.O ) 

901b  F(RmaT  (2513) 

9016  format  (»3»i 

901b  format  <*3  SELECT  EVERY  »,I2,*  SAMPLES*/) 

9019  FORMAT (*  •.F7.?,2x,E12.4,2X,El2,4,2X.F6.?,2x.El2,4,2X.Fo,3,2X, 
1EJ4.4,2X.?6.3,2X.E14,4,4X.F10.4) 

9020  F0»KAT<F13.2.4I5.2F!>.1,  15) 

90  21  FORhAT(//,?X.*W|ND  SPEED* , 3X , .NLC« » 3X , *NAC* , 3 X, «N6 », 3X, «NUHH 1 ST* . 

2 2X.»SHH0NLV.,3X.*HHD1ST.,3X.*IF ILDAT#) 

9022  FORMAT <f9. 2. 1 8 . 3 1 6 , 2F 10 . 1 , ! 10  t 

9023  FoRhAT(lOIR) 

9 0 24  FliRMAT(//,3X.*KFILE*»2X.*LFlLbl*.2X.*LPlLE2*.2K,*LF1tUE3*,2X. 

2*L)  1LE4*.  ?X,»Lf  ILE5*.2X,*LFlLfc6«,2X,*L)'lLE7*.2X,*13HEEH*,3X# 

J*NF|.  A(>.) 

9325  FORMAT (/, 5x,.HRARF(1)**.F7.3.5X,*HBaRF(2)**,F7.3,sx,*TSkECAVQ«*, 

2F3.i),/) 

9026  FuRmAT(/.5X,* FILTERED  DATA  ON  OUTPUT  TAPE*./) 

902/  F0RmAT<2X, ‘RESTART  * • , E 3 . 1 J 

902b  format ( !5) 

9029  F0RMAl(lHi.2X,*NSETS  * » . 1 3 . 5X , *NREPE A T **,I3> 

9060  F0RmAT(///.5X, .APPARENT  FREO» . 5X , * A/C  SPEC* . 5 X . *TR Jt  FWtO*,SX, 

2 ‘RAVE  S»mC«> 

9061  FORmaT(8X,F6.3.7X,E10.4.6X.F6.3,6X.E10.A) 

1925  FORmAT(/#?X,*Y<1R)  • *. 15X , *Z ( IR ) * •> 

1920  FORmAT(iOX.F8.4.15X|F».4> 

192/  F0RmAT(39X,.T1*.8X.*T2*.8X,*T3».8X.*T4*.8X.*T5*) 

1928  foRn at ( /» sx, .frequency  response  functions  ■ «.5nu,s./> 

1929  FORMAT!/. 3X..AVGY  ■* , F 8 . 3, 1 U X . « A VGZ  **.F8.3,//) 

i960  FoRmAT(/. JX,.SPECAVbN(l)  »» , FS  . 2 , SX,  *SPECA  VGM  2 ) *.,Pb,2) 

1963  FoRmAT(3X,.YS».,F10,6.5X.*ZS«*.H0.6)  • 

1964  FoRmAT(2X,.U(1,1)  •• • F 9 . 6 ) 

196b  FiiRmAT(2X,.U(1,2)  *• . F 9 . 6 ) 

1966  FoRmAT<2X,»I  I4,»X .*Ul I ,1 > ** . E10 . 3 » 5X , *U < I . 2 ) **ibiu.3> 

PA3<5E  1 

CALI  PL0T>(PlTARRAV.1U26,1) 

Call  fact3«(i.) 

C ONLY  ONE  DATA  CARD  FOH  NSETS  IS  NEEDED  FOR  EACH  RUN, 

RE*0  9 028,  NSETS 

IF  ANY  NE 2 TaPFS  ARt  NEEDED  WHEN  NHEPE*  f CHANGES  VaL^E,  SET 
1FAUSE  • \ IN  PREVIOUS  FILESXS  DATA  DECK. 

Do  M900  N ?FpE A T * 1. NSETS 

INITIALIZE  PARAMETERS 

IH2mREV  * 0 
HPLOTIlF  s 0 
Pi  = 3 ,14  VS92654 

TK  OP  I * 6,283185308  PSHAH-65 

« * 32 « l’/?5 
GXP 1 * G*  * I 
TrOi  VG  ■ 2.U/G 
GDY1RO  ■ 3/2. 

VSR i P • 0, 

IaSuVEFb) 9mLoRF  «0 
SPEtN*XF=rNUMSPEC«T5PfcCAVG«U. 


BHIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

EBOM  copy  FURKlStiED  TODDC. 
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no  o o n n o o r.  on 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
THOM  COPY  FURNISHED  TO  DDC - 


t; 

C 

c 


Dll  50.  N=t.? 

Du  50.  1*1,4097 

31)  FREuSPEC  l 1 ,N>»0. 

RMj  IN  DATA 


1U0 

2UU 


HfcAi  9031.  NOS.  I M 1 1 > , !M(2)  .UELT.V.CUT.H.NUKPT.LUilHl  ,N,  .NAVGL.UF 
l)  (EOF  ,6U  ) 9999,21)0 

II  * 0 


RfcAo  9032.  LPM1 

Rfc Au  9004.  N1 ,N2, IP AUSf, PH INT.PHAS.PL  rON.NHtCSKP.NOHAyPLT.NPEN 
• .NSPECAV3. IPLOTSP 

HfcAi,  9014,  NREC,  H.LO,PLTSPACE.MTHEH.U»,THfeiA,DKLAStP 

RE  A I 9020 , H.NLC.NAC.NS.NUMHlST.SWHONLY.HHDIST, I)  IlDaT 
RfcAi  9024,  KF  ILE.  UF  ILM  I ),  1=1,7).  ISWtEP.NFLAG 


RE*i  9015.  NGSPEC 
Rf-Ai  9020,  RESTART 


D)Li5t9«  «nr  LASER 
ll«NTOT»NlS 
XmAaVl  * 1.6266 
REX  1 5 ■ 0, 


RaVg  * n. 3 


PRUT  out  CONTROL  CARUS 


215  PH  1 nT  9029.  NSFTS,  NKEPEaT 

PRUT  9004.  N0S,IMm.IM<2>»DfcLT»V.CUT,H.NLMPT.L3WHl,NF,NAVGL» 

• D) .t  9H1.NI.N2,  IPAUSt, PRINT. PHAS.PL TON, NHECSKP.NOHAVPLTiNPtN.NSPtCA 
•Kb. i PLOTS*, NREC,  LI.LO, PLTSPACt , MTHEH. UP, THE Ta.DK.aSEH 
PRINT  9021 

PRUT  9022.  W.NLC»NAC,N3»NUHHIST,SWHONLY,WHClST.lri)H)Ai 
PRINT  9024 

PRINT  9025.  KFILE.ILT ILE<J),)*l,/>,!SWfcEP,NFLAG 
PRINT  9015.  NGSPEC 
PRUT  9027,  RESTART 
I f < i N( 2 ) , N6 , lM2PRtV)  220. 26U 
220  IFS  ■ -1 

I M2PREV  = I M ( 2 ) 

260  IfUPEN  .ji).  0)  2/0. 2/5 

2/0  NpEn  s i 

DOHNMUD  TRACKS  REQUIRE  NEGATIVE  (•)  AIRCRAFT  VELOCITIES 
CONVERT  AIRCRAFT  VELOC11Y  FROM  KNOTS  TO  FEET/SEC. 

If  OMEGA  ■ theta,  DUPRLPR  SHIfT  SPECTRA  AS  If  MAVtS  MuVbD  PARALLEL 
TO  WIND.  IF  OMEGA  ■ U OR  ISO,  IT  IS  ASSUMED  THAT  4aVES  MOVED 

parallel  to  flight  track. 

2/5  OMEGA  ■ T4ET4 

If (THETA. iT. 90. AND, THETA.LT ,2/0)  276.27/ 

2/0  OMEGA  * 13U. 

GO  50  278 
2/7  OMEGA  « 0, 

2/8  V * V«1.69/8*C0SF(PI»0MEGA/18u. > 

If ITHETA. jO.9o.OR.THETA.EO.27U)  280.290 
280  V«0. 

SKIP  NRECJKp  FILES  ON  OUTPUT  TAPE  OR  ON  INPUT  TAPE  |F  Li  ■ LO. 

2VU  DO  -'50  I « l.NRECSKP 
CALL  SK1P--ILF  (LO) 

55U  CONTINUE 
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ORonn  non  noon  oonnci 


IM  VSHIP.  =0, 0.  ) 5 2 3 » 5 w;  0 


CALCULATION  WHICH  ALLOWS  SHIPBOARD  hadak  and  aihbotne  phqulometer 
T L ANAL  Y7E  DATA  OVER  TRACKS  0^  SIHlLAH  LENGTHS.  IT  UtiERMINtS 
The  PRO»>-R  sample  interval  and  EVERY  ISEC-TH  POINT  WILL  BE  USED. 

520  TOShTR  a inZ4,*VSHIP*U,125 
VA0S  * A85f<V) 

TlRrAF  a TflSHTR/VABE 
RlHl  a IM(i) 

TjRhDA  * TIREAF/RIMl 
RlStC  a TIWEDA/DELT  ♦ 0.99 
IV  ( RISE:  .LT.  1,  ) 521,522 

521  RlSi-C  a 1, 

522  I SEC  a RIS^C 
DELT  a IS?C*I)ELT 
GO  13  5?4 

IT  0 DESIT’-D  SAMPLE  INTERVAL  IS  LARGER  THaN  THAT  OF  THE  ORIGINAL 
DATA  set  or  than  lHAT  DETERMINED  ASOVt,  cetine  the  UtSlREU  UELT 
IN  the  INPUT  DECK  AND  THE  BELUW  EQUATION  WILL  PRDPEHLY  DEFINE 
I SEC  PRJVIDED  THE  DENOMINATOR  is  correct. 

523  l Sfet  a Dfc.T/0.nilX3J3 

524  PR  I « T 90T  5 ■ ISEC 
PSEl  s I S = C 

PLTsPACE  • PSECaPLTSPACE 
OX  a ASSf  t V.OElT ) 

NUMRTPl  * NUMPT  ♦ 1 

Tt-SI  IF  FILTER  IS  TO  BE  COMPUTED,  IF(NUMPT  .GT.  U>  » YeS 
52b  IT (NJMPT.sO.O)  1800.550 
COMPUTE  f li.TER 

NOMpT  * NJHBER  OF  POINTS  OVER  WHICH  THE  AVERAGE  IS  TO  Bfc  TAKEN 
CUT  a CUT-OFF  FREOUENCY/SAMPL 1 NU  FREQUENCY 
H a S .OPE  OF  WEIGHTS 

550  KA  a NUMPT  ♦ 1 
Cl  M a CUT  * H 
FILI(KA)  1 2 . *CTM 

SUM.  a 0. 

DO  /00  I a 1, NUMPT 
P ■ I 

QO  * 1,  - (14.«H*H»PaP> 

IT (UQ  ,Nfc,  0. ) 600,57b 

5/5  FILT(I)  = it, 

GO  10  700 

Su  0 FILM!)  a (COSF  ( TWOPI*P*H)*SINF(  TWOPUPaCTH)  )/(Pl  *s*QO) 

SUMA  a S04K  * FILTH) 

700  CONI INUE 

FLJl)  « (1,  - (FILT(KA)  ♦ 2 , *SUMK  1 ) / ( 2 ,*  NUMPT  ♦ 1.) 

DO  p 00  I a 1 , KA 

WfclNTl  I ) a F ILTI  I ) ♦ TLJO 
So 0 CONTINUE 

If IlDWHI . -0,2)  825, B50 
325  WEl(.T(KA)»  1.  - WEIUTtKA) 

DO  *40  lal, N'JMpT 
Sail  NElGTUja  -WE  I GT  ( I ) 

C CAtlULATc  A Nil  PRINT  FREQUENCY  RESPONSE  FUNCTIUNS  0?  THE  FILTER 
050  DO  *00  H»ai,S 
BR2aFwFACrOR(MP) 


IBIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
ISOM  COPY  FURNISHES)  TO  DDC  - — 
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orsoortr:  an  an  a 


? 


9UU 


lSOU 


FR*| ,/(|H(ll *BP2> 

FORi  SF  ( RP ) *mE 1 GT  < KA ) 

DU  '.00  UJ*1,NUMPT 
FOR.  SFIrP)*FQRFSF(MP> 
PRMT  192  7 

PhliT  1923. (FQRESFIMP) 
CUN  i IVJfc 


BUS  PACE 

kbukoof* 


“ sHssioms 


♦ 2.«WE|GT( |KJ>*CUSF<TkUPl*l<J»FRj 
. HP»1.5> 


R t-  A > IN  U 4 T A 


IMl  I ,E0,  LO)  1810.1890 

if  wcadinj  in  more  than  one  file  per  run.  no  need  to  remind  u. 

MF  Hf  LY  ADJUST  NRECSKP  ACCORDINGLY  IN  D*  I A DECK  OF  *ULL0mING  FILE. 
1810  BuF> EH  |N(l(.1)(Y(1),Z(6192)) 

1815  IF(uNlT.Ll)  1815. 1820. 1820. 1820 
1920  Call  SK|P?4LE(!  I ) 

GO  13  219 J 


UnP.CK  and  assign  values  to  y and  z 

for  un»ac<ing  files  in  seuuence  off  a tape,  cu  not  rewind  li. 

MERi-LY  SET  KF|i.E*l*NUMBl:R  Of  LAST  FILE  UNPACKED  IN  U*T*  DECK  OF 
NEX1  FJlE,  BUT. IF  LI  IS  REWOUND,  KFlLE  MUST  * 1 INSTEAD, 

19*0  II  = 0 
NP*1 

DU  >155  I f I lE  a 1,'TF 
DO  1870  JZ  *NP , 7 
IUlFILEI  J7>  .EO.O)  18/0,1860 
I960  MF 1 l E*  JZ 
Nw» jZ*l 
GO  10  1883 
19/0  CONTINUE 
i960  call  rota»e 

DO  2150  MFC*  l.NREC 
DO  2150  I SMP  = 1, ISMEEP 
CALi  UNpC<(  1SWP,  IREO.NFI.AG) 

11  * 11  ♦ 1 

Y ( I 1 ) ■ :HAN<NLC>*UFLASEH 
Z(ll)  « CHAN(NAC)*DF 
2130  CoNl INUE 

Call  s<|R?!lf<l  I ) 

2135  CONTINUE 

C OPTION  TO  PRASE  SHIFT  DATA 

IF (PMAS.EJ.O)  2165.2160 
2160  CALL  PHASr(NTOT.QP) 

C OPTION  TO  RECORD  DATA  BEFORE  UR  AFTER  FILTER 

2165  IF  ( IFILOAT.EO.I  . A NO , NUHPT . NE . 0 > 2190.21/0 
21/0  BUF.ER  OUT  < L0 . 1 > ( Y < 1 ) . Z ( 8192 ) ) 

21/5  If IuNIT.lD)  2175.2160 
2160  END!  ILE  L 3 

IF  ( IFILDAT.EO.1  .AND.NUMPT,  NE.il)  2185.2190 
2165  PrInT  9025 
GO  TO  399J 

2190  IF  (mTREM.  3T.o)  2195. 220(1 

2195  RERINJ  Ll 

2200  NIDI  * NO 5 * II 

C 

c print  f,rst  value  of  each  record  uf  v and  z. 
c 

PRINT  1925 
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n n n n non  non 


DU  2205,  12  = l.NREL 
Ih  = (12  • l)*ISWfeEK  ♦ 1 
2205  PRImT  1925.  Y(  IRt.ZlIH) 

Cali jl  atc  ym4X.ymin»zmax,  and  zmin. 


2212 


2220 

2230 

22*0 

2250 

2260 

22/0 
2260 
2 JO  0 


23/6 


YmAx  * JMAX 
YkIn  * /M|\ 

pr^t  930  s. 

Du  -3)0  I 

IV  | V ( l 1 . JT 

YMAX  s Y<  I > 
GO  10  226 J 
IKrm  . .T 
YMlN  a Y<  I ) 
IF  (ill)  • 3 T 
ZMAx  3 M) 
GO  10  2300 
1M/U)  ,.T 
ZHi  s Z(  I) 
CONI  1 NUE 
PHUT  poo;. 
PRl\T  9007. 
PHUT  9014 
U (mExTS. ;o 

in  i stc.Gr. 


* -100, 

3 100. 
nTOT 

3 1 , NTOT 

, YMAX)  2220. 223U 


. YM I N ) 2240. 7250 

, ZMAX)  2260,2270 

.f.  ZMIN)  2260, 2300 


YMAX, TM1N 
7MAX , ZM IN 

,0.)  2376.341)0 
1)  ?300<  >*440 


StLi-CT  LA$tR  AND  ACCELEROMETER  DATA 


2360  II*) 

DO  >600  1 = 1 ,NTQT  , 1 SEC 
Yl 1 ) • Y(  1 ) 

2(1.)  a 2(1) 

2600  II  = II  *1 

JIS*-C  a I l ♦ 1 
I) ( iSfcC.EJ.i)  2660.2610 
2610  Du  >650  1 a JISEC.NTOr 
Y ( I > a 0 . 

2 650  £ 1 1 > 3 3 . 

2460  NlU'aNOSaJI 

fin  average  values  before  filtering 

Go  10  344) 

AMP. Y FILTER  TC  data  (0  a NU) 


2700 

2400 


2900 

3100 


RaVi.  s 1.) 

N1 0 1 3 nTJT  - ?*NUMMT 

MIN  a 1 

Do  3200  I 3 MIN.NToT 

XH  s I * NUMPT 

SUN.K  a WMGT(KA)*Y(KH) 

Sum.  < a hMGT(ka)*Z(RB) 
do  >900  J a l.NUMMT 
Ia  s NUMPT  ♦ 1 - J 
IR  s U ♦ NUMPT  ♦ I 
It  = J - t ♦ I 

SUMCK  a SJMCK  * UE|UT(lA)a(Y(lC)  ♦ Y<IB)> 
SUMl.X  a SjMUK  ♦ Wfc  I UT  ( I A ) * ( Z ( I C ) ♦ Z t I B > > 


CONI  I NUE 
Y( I ) a SU4CK 


practxcaw* 
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raou  cqpx  FyjpistiEP  °^c-  - — 


Z(  I ) = SU4I  K 
320  0 CONTINUE 

C SET  Rfc  S I •)  j a L DAT*  PTS  To  0. 
inorPi*NTjT  ♦ i 
DO  .'250  1 «NT0TP1,N06 
325U  ?( l )*Y(  1 ) > ' . 

C RECOMPUTE  t-XTHEmES  UF  T AND  7 AT  TER  FILTERING 
RfcX  l S * i , 

GO  13  2212 

34uu  IM  |F|LD»T.E.a.l)  21/0*3440 
C 

C F i N AVERAGE  VALUE  UF  INPUT 

c 

3440  Avar  * A V 37  s 0, 

1TMAVGL  .NT.  NTOT)  3450,34/5 
3430  NaV  ,l  * NTOT 
34/5  Dfj  '5UO  . 3 l.NAVbL 
AVGV  s A V j Y ♦ Y(I) 

AVU2  s AVJ/  ♦ 7(1) 

3300  CONTINUE 

AVGy  s a'/JY/naVGL 
AVU7  * AVJ7/NAVGL 
PH  I \'T  1V29,  AVGY, AVGZ 
IT  (RAVG.fc J. U . ) 3550,3575 
3330  IM \UMPT . iO. 0 > 3575,2/30 
35/5  Do  .'600  I s l.NTOT 
2(1)  = 2(J)  - AVGZ 
35u0  Y(I)  s Yd)  - A VGY 

YHAx*LABU/E(1 TrYHAX  - AVGY 
YhU«l.8fc!.DW(l>=YMIN  - AVGY 

ZrAX3LABUVM2)3ZMAX  - AVGZ 
ZMI'«LREl3W(3)=ZH1N  - AVGZ 
NSPtCTRA  « NTOT/IM<l) 

C 

C IT  (PlTQN  .GT.  0)  UO  NUT  COMPUTE  SPECTHA 

C 

IT  (►  LTON.  3 T . 0 ) 3625.3050 
3625  JlOlsNToT 

NbTMSTsNTxNTlMfcS»l 
NbP*  CAVG  « 1 
GO  TO  43/3 
3630  II  = 1 

12  3 I H ( 1 ) 

1 4 3 12/2  ♦ 1 

1 6 3 I 4 - 1 

OECPtS  = 1 ./<7,*DELT*IM<1> ) 

FaCo  3 pt.T*!M(l)/TNOPI 
Nil.-  ES«MS»FCTRA/NSPtCAVU 

C FOR  THE  UJLF  STREAM  EXPERIMENT.  SET  NTJMES  * 1. 

NT  I "ES  3 1 

It ('TIMES  .LT.  1)  3660,3675 

3660  NT  1 ciES  ■ 1 

C IT  iv  T 3 'i  T I 4 E S CORRESPONDS  TO  DATA  HH  l CH  SHOILC  NOT  3E  ANALYZED  AN* 

C It  Nlal , SET  NTsNTIMES  - 1, 

36/5  DO  t 450  NT  « 1 , NT  I MES 
Do  4 150  N 1 1.,? 

SpE  T A VO  s 0 . 

I A Dl  3 p 
NUTT  s u 

DO  '700,  1*1,5500 
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J7u0  U l 1 , N > ■ 3 . 

NSThIST  * (NT-1  )*NSHEI;aVQ*!H(1)  ♦ 1 
JlOl  S NSTH1ST 
NS  Tm  = 0 

Du  -ion  i * i.nspecavg 

C Tk-  S I r OR  3001)  DATA.  NSTH  I ST»S1'ANT  OK  DESIRABLE  DA  I * AND  LOCATION 

C HmERE  HISTnQRAM  BEGINS.  FOLLUWlNU  DO  LOOK  ALLOWS  OS  lO  KICK  U'.'T 

C St-  SMENT  j OF  LENGTH  = 1 M ( j. ) 1 RUM  THE  ANALYSIS. 

Du  3730  J * t.15 

If ((NT-1>»NSPECAVU  ♦ 1 , EQ , NbSPEC  ( J ) ) 3705,3720 
3 7 U t>  NSTh  =N5l  4 * 1 

I) (NSTH.E3. I ) 3710,0715 
3 7l U NSTmIST  3 itlM(l)  * NSTH 1ST 

3 7 lb  Go  10  410  3 

3730  CUNlIN'Jfc 

N10T  s NTJT  * I M ( 1 ) 

L 3 1 

1a  = ( (Nl -1 )«NSPECAVG  ♦ J . 1 ) * 1 H ( 1 ) * 1 
lb  3 ( (Nl -1 >*NSPECAVG  ♦ 1 ) * I R ( 1 ) 

I) ( ( Id-IAJ .EO.NAVGL)  3758,3751 
C F 1 Nn  AVERAGE  OF  THE  DATA  SEGMENT  (IA.IB) 

3731  YS  3 ZS  3 11. 

IMN.EQ.l)  3752,375- 

3752  DO  ,-753  J Y 3 I A , IB 

3753  YS*TS  ♦ Y(.)N  ♦ 1 ADD ) 

YS  3 YS/12 

GU  10  3755 

3734  DO  .1755  JY3IA, IB 
375b  ZS»/S  * Z ( JN  * I ADD ) 

Zb  3 ZS/I 2 

3756  PRINT  1933,  YS.ZS 
373b  DU  3800  J 3 IA.IB 

If  ( N .Eli.  1 ) 3760,3770 

3760  0 ( L } s VC  J * I A D n ) - YS 
Ud  70  3003 

37/0  OIL)  = Z C J * I A D 0 ) - ZS 
3 8 11 0 L * L * 1 

call  FOURIER  <G),S.  Ih(2>.  If  S) 

38/0  RoS  ■ 2 

C CALCULATE  POWER  spectra 
3BV5  DO  3900  J 3 1,14 
39u0  O(J)  3 0(  2-.J-1  > **2  * U(2-J)*«2 
DO  «000  J«l, 14 
4000  U(J,N)»U<  J,N)  * 0< J) 

SPE.TAVG  • SPECTAVG  ♦ 2.0 
If (N.bO.l)  4060,4070 
4060  PRINT  1934,  0(1) 

GO  10  4103 

40/0  PRINT  1935,  0(1 ) 

4100  CONI INUE 

SPfcl AVGN( Y > 3 SPECTAVU 
4130  CONTINUE 

TSPt-CAVG  t TSPECAVG  ♦ SpECA VGN ( 1 ) /2 . 

C 

C TEST  FOR  END  OF  GOOD  DATA  SO  AS  TU  AVOID  PLOTTING  JNDES 1 RED  DATA 

C AT  THE  E Yl)  OF  DATA  RbCURU  HEiNU  ANALYZED  IN  ABOVE  LOOPS, 

C 

Nf  IN  3 NT aNSPECaVG 
JSTrjP  a 0 

If (NGSPEC(MFIN).eQ.O)  4154,4151 


oc:o  noon  n n c: 


OJIS  PAGE  IS  BEST  QUALITY  FRACTICABljB 
*58QM  OUPY  P^RJUSHED  TO  HDQ  - 


4191  NS  T A a (Nr-1)*NSPECAVU  * i 
JSTOP  a 1 

DO  4163  HP  * NSTa.NFIN 
JJP  a NF  H - IMP  * NSTA 
JjPMl  a JJP  - 1 

lF(NGSPfcC< JJP) .EQ.NUSPECI JJPM1)  * 1)  4162, 4l54 
4195!  JSTOP  a JSTOP  ♦ 1 

4193  CONTINUE 

4194  JlOIa  JTOr  - NSTHISY  * ( NSPECA VU- JSTOP ) * I M 1 1 > 

PRINT  J.93D,  SPFCAVGN(1),SPECAVGN<2> 

IF  (SPfcCTArfG.EQ.O.)  641)0,4156 

APPLY  HAM4ING  TO  THfc  TOTAL  VARIANCE 

4196  DO  4170  Nal,? 

Ud.NJa  U,54*U(1.N>  ♦ 0.46*0(2. N) 

DU  *160  i *2,  16 

4160  U(l.N)a  0,?3*<U(I-1  »N)  * U(I*1.N>)  ♦ 0.54*U(1.N> 

41/0  U«l4.N)a  D.64*U(14.N>  * 0.46*0(16. N) 

uo  loop  to  transform  acceleration  to  displacement  spedtrum 

IRANbfORMATl DN  SHOULD  ONLY  BE  INITIATED  FOR  ACCELERATION  DaTA 

41/1  U(l.l)*u( t,2>»0. 

DO  41/5  L4*2,  14 
Re2*ILW-l ) 

FTRUaR.FHrORES 
F 1 R u 4 a ( .'ITRU*TW0PI  >**4 
41  /i>  U(Lh,2)eU<LW,2  VFTRU4 

AVERAGE  T4F  VARIANCES  AND  CALCULATE  SIGNIFICANT  NAVE  HElQHT 

EMAX(l)  * PMAX(2)  ■ 0. 

41BU  DO  4200  1 a 3,14 

U ( 1 , 2 ) * J ( I , 2 ) /SPECA VGN ( 2 ) 

U(I.1>  a J( 1 ,1)/SPECAVGN(1) 

C N3*l.  SUBTRACT  A/C  MOTION  SPECTRUM  FROM  LaSE*  SPECTHUM, 

IF  U3.EQ.  1>  41B2.41V0 

C TfcSI  VALUES  OF  THE  VAR  I ANCES-- DU  NOT  ALLOW  NEGATIVE  VALUES  OF  U(1> 

4162  IF  (U( l ,1> ,LT.U<  1 .2) ) 4164,4106 
4164  PRUT  1934,  I . II  ( 1 , 1 ) . U ( 1 , 2 ) 

U< I ,1)»0  • 0 
GO  10  4193 

4166  U( I , 1)  a J(l,l)  - U ( 1 , 2 > 

4190  EMAx(l)  a EMAX(l)  ♦ UU.l) 

EM  A X ( 2 ) a R M A X ( 2 ) ♦ U( 1,2) 

4200  CONI INUE 

4203  HHAR(l)  a 4,0*SQRTP (EMAX(l) ) 

HbAR < 2 ) 1 SQRTF (EM'X(2) ) 

C 

C WRITE  HEADINGS 

4210  PRIM  9003,  LPHl.NTOT.OFLASER, IM(1),CUT.DELT,H,V.NJMPT.EMAX(1), 

• HBl  R ( i > , iMAX ( ? ) , H« AR 1 2 ) . 1 F 1 LT ( LOWH l ) 

4?16  HBARF<l)a4RARF(l)  * HBAR(l) 

HbARF ( 2 ) a HRARF ( 2 ) ♦ HRAR( 2 ) 

TnUpSPEC  » TNUMSPEC  ♦ 1..  0 
IF (M .FQ. I.AND.NT.EU.NTIMES)  4216,4217 
4216  HbARF(l)aHPARF(l)/TNUMSPEC 
HhANF(2)84BARF(2)/TNUMSpEC 
421/  IF (SWMONLY.EO.l.O)  4325,421V 


C UO  LOOP  FOR  JALCUIATION  OF  POWER  SRECTHA  AND  RELATED  3U*NTiTjES 

c a DOPPLER  shift  algorithm  is  included  to  remove  the  fReouency  shift 

C DUE  TO  THE  HORIZONTAL  AIRCRAFT  MOTION.  IF  N?  EQUALS  1.  aN  AUDITION 

C aL  Cl'RHECT I ON  IS  APPLItD  TO  ACCOUNT  FOR  ANGULAR  SPHEaDiNu  AS  PRE- 

C DlcTcD  8Y  si-ILLIP.S  RESONANCE  THEORY  AND  JS  APPLICABLE  F OR  FLIGHTS 

c imat  paral.el  the  hind  direction. 
c 

42iu  on  -3un  n ■ 1,2 

II (N.feO.2)  4280,4220 
42211  RRlr.T  90  0 A 

4221  MSA  I n i 

C HaNT  10  PRINT  every  nz-th  pi  PLUS  FIRST  So  RT, 

N2  * 14/64 
ANZ  * NZ 

DO  «265  1*1.14 

R 3 2«  < I * 1 ) 

MRU  = R*  •■REORES 
IKY.EQ.C.)  4222,4223 

4222  FmuhT»2.*»I*FTrU 
GO  TO  4224 

4223  CONST  » G/(2.*V) 

CuNSQ  * LDNST*C0NST 

FMUHT  ■ SJRTF (CONSQ  ♦ GXP l *R/ ( IM< 1 > *0X > > - CCNST 

4224  IMl.EQ.l,  AND.  CoNiT.  (iT.  o.)  4225.4226 

4225  FhUhT  * 0, 

4226  AS 01  * U l I , N ) *DELT* 1 M ( 1 ) 

IMV.EQ.O,)  4228. 424(1 

4228  ASO^Ol  » ASQt/TWOpl 
GO  TO  4231 

4230  ASQMU1  * AHSF(TU0DVG*(FMUHT  * CUNST)*aSU1«V)/ TWOP1 

4231  AS0A  T 1 ■ ). 

1MN2.60.D  4232.4242 

C STATEMENTS  4232  TO  4242  PERFORM  THE  RESONANCE  THEORY  DOrPLER  SHIFT 

4232  IF (I  HURT. =0. 0 . ) 4234,4236 

4234  F*T  • ASO  403  * ASQ2NIK  I ,N)  * FNUWT  ( I «N.)  » 0. 

GO  TO  4244 

4236  FRUf-T  ■ H*FMUHT*«2/U 

4240  ASQHU1  ■ AS0MU1*SQRT(U/(FMIIHT*R)  J/2. 

4242  FNUhT(J.N)  a FMUHT/TWOPl 
TKT  • <FMJHT**2»/G 
AS02NU ( 1 > Y > * THOP I • ASQMU1 
FREoSpECU.N)  = FREUSDECI  1,N>  * AS02NUU.Nl 
IF <F HURT. NE.O , ) 4244,4246 

4244  ASOxTl  ■ ASUMU3 •GDVTWU/F  MUHT 

4246  If  (PRINT.  50. 0 .)  424'/, 4265 

4247  IF < I4.GT.l28)  4248,4249 

4248  RRNUM1  * I 

PRNUN2  ■ »RNUMJ /ANZ 

1RRJ.  * PRNUM2 

PRNUM3  ■ IPR1 

PHNUH4  * >RNUM2  - PRNUM3 

IF (I .LE.50.OR.PRNUM4.EO. 0. ) 4249,4265 

4249  IF ( | .LE.NJMPT  .♦  1)  4290,4255 

4250  PR  1 nT  901».rTRU,U( J # N) , ASQ1 , FMUHT, ASQMOl . FNUR T ( I . N ) , ASQ2MUJ I , N 1 . 
•FKT , ASOKT 1,HE IGT( I ) 

GO  TO  4265 

4255  PRINT  901>.FTRu.U< IfNI.ASOI, FMUHT. AS0MU1,FNUFT( I,N),ASB2NUU.N), 
•FkT.ASOKU 

4265  continue 
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CO  10  430) 

4200  PRINT  901) 

DU  4290  I a 1,14 
R * 2* ( l - 1 ) 
riRu  a H**PE0RES 
FNUhT(  I,N)  a FTRU 
FMUmT  a f fRU«TW0Pl 
ASO,  « U<  1 >N)*DELT*IM<  1 ) 

ASU2NIJ2  a AS02 

A&0MJ2  a ASQ2/TW0PI 

If (PRINT. sn.O. ) 4280, 42b9 

4205  IfU4.GT.i28)  4286,4288 

4206  PHNUM1  • I 

PRNUM2  * »RNUH1/ANZ 

IPRl  s pR VUM2 

PRNUM3  a JPR1 

PRNUM4  a »RNUM2  - PHNDH3 

If ( 1 .LE.&).OR.PRNUM4.fc0.o. > 4280,4209 

4208  PRIM  90U,  FTRU.UU»N>.ASQ2.f  HUMT.ASQMU2 

4209  ASQ2NUC 1,4)  * ASQ2NU2 

FREOSPEC  < I • N ) = FREUSPECI  1 ,N)  ♦ ASQ2NUU.N) 

4290  CONI INUE 
4300  Continue 

calculate  rave  height  distribution  function 
If  SPECTRA  (1PL0TSP  NE  U)  ARE  NUT  Plotted,  teen  PLOT  fOH 
THE  WAVE  HEIGHT  0IS1RI8UTI0N  MAY  BE  POSITIONED  WRONG,  ALSO. 
ThEhE  IS  VO  EQUIVALfcNT  10  FREUSPEC  FOH  I M I S FARAMtTbR. 

4325  HO  6000  N * 3,2 

I aUi;  s o 

43/S  Dll  4400  M a 1,301 
4400  DfNi(H.N)  = 0, 

Du  4500  M a 1 , NSPECA VG 

Ha  = ( (NT-! )«NSPECAVG  ♦ H - 1)*|M(1)  * 1 

MB  = ( (NT.n.NGpECAVG  ♦ H>*IMll> 

DU  4415  «M  a x , 25 

II  <(NT-i J.NSPECAVG  * M .EQ.  NUSPEC(MM))  4500,4415 

4415  CONI INUE 

UU  *•5110  MM  a HA, MB 

If  ( N .=0.  1 ) 4416,441V 

4416  INDcX  a Y < mm* I ADD ) *4 . 0 * 151. S 
GU  10  4415 

4417  INDcX  a Z ( mm* 1 ADO ) »4 . 0 ♦ 151.5 
4410  If  UNDEX  ,LT . n 4425,4450 
442b  IljOrX  a 1 

GO  TO  4485 

4450  IMiNUEX  , GT . 3(H)  44/5,4485 

44/5  lMli-X  a 331 

4405  Df Net  INDEX, N>  a DFNC ( I NUEX , N ) ♦ 1 
4500  CuNi INUE 

DO  a 60  0 m z 5,301 

IfCuFNC(M.N)  , GT . 0,1  4700.4600 

4600  CONI INUE 

4700  IST(N)  a i-XMOUF(M,4i  ♦ 1 

If ( I ST ( N ) .LT.  1)  4/25,4750 

4725  IbT(N)  a l 
4/50  Du  4800  MH  a 4,301 
(1  ■ 301  - MM 

If  (l>FNC<M,N)  .GT,  0)  4900,4800 
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4800  CONI  IN'JE 

4900  INU(N>  a 4 ♦ 1 * XMOUF ( M , 4 ) 

II  ( 1 ND( N T .GT,  301)  4V25 , 4950 

490b  ] NO « N > a 301 

4950  CENTROID  • STANDDEV  a SKEWNESS  a KUHTOSlS  a 0. 

DtV  a ( I ST ( N ) - 151I/4. 

I » IST(N) 

J « IND(N) 

FNTOT  a N TOT 
DO  boon  M a I,J 
CENTROID=:ENTROID*0lNC(M,N>*DEV 
SlAf'DUEV  a STANDDEV  ♦ Dl  NC  (M,  N ) *DEV**2 
SKEwNeSS  a SKEWNESS  a Dl  NC  ( M , N > *DE  Va*3 
KURTOSIS  a KURTOSIS  ♦ Dl NC(M,N>*0EVa*4 
bOOO  DfcV  a DEV  * 0.25 

CENTROID  • CENTROID/FNTUT 
SIANDOEV  a STAND DEV/FNTUT 
SKEWNESS  a SKEWNESS/TNTOT 
KURTOSIS  a KURTOSIS/FNTOT 

KURTOSIS  a KURTOSIS  - 4 . *CENTHO( D*SKEWNESS  ♦ 6 . *CE NTRO I U»CENTHOl na 
* STANDDEV  - 3 • ‘CENTRO l D*« 4 

SKEWNESS  a SKEWNESS  - 3. ‘CENTROID* STANDDEV  ♦ 2 .aCE NT«Ol U*«3 
VARIANCE  a STANDDEV  - CENTROlU*CENTROlD 
SIAnDDEV  a SORTFIVAKIANCE) 

SKEWNESS  a SKEWNESS/(STANDDEV*VAR1ANCE) 

KURTOSIS  • KURTOSlS/<VAHIANCE«VARIANCE> 

PRINT  901;.  N.I.NiJ, CENTROID. STANDDEV. SKEWKESS.KJHTOSjS 
600(1  CuNl  INUE 

MIN  a XM1  NOM  I ST ( 1 ) il ST < 2 > > 

MAX  a XMAxiiF(|ND(1)iIND(2>) 

IbElOW  a (154  - MINI/4 
I ABOVE  a (MAX  - 1481/4 
IbElOW  a IBELOW  - 4aM0D( IBELOW, 2> 

1 ABOVE  a I ABOVE  * 4*MUD < I ABO VE , 2 > 

MIN  a 151  - 4a  I BELOW 
Max  a 151  ♦ , 4a  I ABOVE 

IMmIN  . L 3 .AND,  MAX  ,GE.  2”>  6025. 6030 

6025  IbELOW  a JAfinVE  a 50 
6030  II (HBAR(l).LT.O.)  6070,6090 
6090  II  ( 1 BEL0W.*i.LT.  IBELOW)  6|i55.6060 
605b  IbELOWF  a IBELOW 
6060  II  ( lABOVE.'ULT.  IABOVE)  606b, 6070 
606b  IABOVEF  a t ABOVE 
C 

C PLOT  DESIRFD  RESULTS. 

C 

60/0  II (NOWAVP.T.EO.l)  6150,6100 

6100  CALL  PLOT  2 ( JTOT,NUMPT, IMID, IPLUTSP.LBELOW.LAbOVE, NUMHIST.N1, 

2 NSTHIST. nT.NTIMES.NOS.SWHONLT) 

6190  II (pLTON.NE.n.OR.SWHONLY.EQ.l.O)  6500,6160 
6160  FSPECMAX  a FDMAX  a 0. 

DO  6400  N a 1,2 
6164  DO  6200  1*1,  14 

II (Nl.eO.l.AND.NT.EU.NTlMES)  6165.6170 
616b  FREqSPEC  ( I » N )aFREOSPEC.(  1 , N ) /TNUMSPEC 
61/0  II(aSQ2NU( I ,N) .GT. FSPECMAX)  6175,6200 
61 /b  FSPECMAX  * A S02NU ( 1 1 N ) 

6200  CONTINUE 

II (SPECMAXF.LT. FSPECMAX)  6215,622b 
62lb  SpECMAXFa"6PECMAX 

THIS  PAGE  IS  BEST  QUAIitTY  PRACTICABLE 
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622b  DO  6300.  1«MIN,MAX 

IfUJFNC(I.N)  .GT.  FUMAX)  6250 . 6300 
6260  FDHAX  ■ D-'NCM.N) 

6300  CONTINUE 
6400  CONTINUE 

IMM.EQ.1.AND.NT.E0.NT1HES)  6403.6406 
C ASSIGN  VA.UES  TO  FNUHTF 
6403  DO  6403  N3«l,2 

DO  6403  I 9»1, 4097 

6405  FNUhTF(ID.NB).  FNUHTUD.NB) 

6406  IMMBAR(l)  .LT.O.  ) 640V, 6407 

6407  111  lPL0TS»‘.EQ.0>  6408.6409 

6406  CALL  PLtS»EC  (AS02NU,f nuht.lphi.hbah, m.i.fspecnax.nohavplt, 
•XMAXVL, ISjC. IBELOW, I ABOVE. MSAT.N2, THETA, Nl.NCS.PLTSPACE! 

640V  If (Nl.EQ.t.AND.NT.EU.NTlHES)  6410,6415 
6410  PRINT  902}'.  HBARF(1J,HBARF(2)«  TSPECA VG 
PRINT  9033 
DO  6412  <Mml,75 

6412  PRINT  9031.  FNUHTF  I KH , 2 1 , FHEOSPBC  < KM,  2 ) ,FNUf  TF  < XM,  \ ) , 

2 FrEQSPE:<KM,1) 

call  plts»ec cfreospsc, fnuhtp. lphi.hbarf. h.i.specmaxf.nuwavplt, 

2XMAXVL.IS5C, I BELOWF, I ABOVEF , MSAT , N2, THKT A.Nl, NOS. PLfSPAUE » 

6415  IM 1SEC.E3.1.AND.HHUIST.EQ.0.U)  6425.6450 

6425  Call  PLOT3F  (FDHAX. IBHLOM, IABOVE.MIN. MAX. LPHl.MPLOTUf) 

6450  CONTINUE 

6300  If ( I PAUSE  .NE.  0)  6600,6720 

6600  pause 

6720  If (RESTART. to.  1.  ) 6900,100 

69uo  continue 

9999  CALL  STOP»LOT 
END 

SUBROUTINE  PLOT2  (NPS.NUMPT, IM. PLTONL Y. LBELOE . LABU YE, NUMHI ST, N1 . 
2 NST, NT, NTIMES, NSAMP, SHHONLY) 

DIMENSION  LBEL0R(2),LA8(IVE(2) 

CqMmON/SF/NPEN.PLTSPACE.DELY 
C0Mh0N/A/T<8l92) ,Z(«192) 

TYPE  HEAL  LABoVE.LBELoR 
TYPE  I NTfc  3ER  PLTONLY, 1 DEC 

C NOTE  that  belt  and  PLTSPACE  HAVE  SWITCHED  LOCATIONS  AND  VALUES, 
210  FCiRmAT(/,  JX,«NSAMP  ■ *,15> 

PRINT  210,  NSAMP 

C FOR  GULTi  STREAM  EXPERIMENT,  NUMPT  ■ 5l2, 

NUMPT  ■ 51? 

Nf Inal  » nst  ♦ nps  • 1 

XXFmAX  ■ NSAHP*DElT 
NPTS  » IM 

Call  nehP:N<npen» 

If (NSAMP  - 2«NUMPT,fcQ.NPS>  3,1 
C RECALCULATE  IABOVE,LBELOW 

1 LHElOW ( l > • LBEL0H(2)»0. 

LaBoVE(1)»LABOVE(2)"0. 

DO  200  I * NST . NF I NAL 
If (LABOVE(l) .LT, Y( I ))  160,184 
160  L4B0VE(1)*Y(  I ) 

Go  TO  168 

164  If  (lBELOWU).QT.V(IM  186,188 
186  LBEi  OM(l)«Y(  I) 

188  If <IAB0VE(2>.LT.Z<!>)  190,192 
190  LABuVE(2)«Z(  I) 

GO  Tq  ?00 
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192  IMLBEL0H(?).GT.Z(J))194,20U 
194  LBEl  OW ( 2 > «Z ( 1 ) 

200  CONTINUE 

3 If (PLTONLY.GT.O. or. swmonl y • to • 1 1 0 ) s,4 

4 CALL  PLOT(0.0,  15.0.-3) 

8 CALL  SYMBOL  ( -1 . 0 , 1 . 70 , U . 35, 19HMEASUREMENTS  IN  HKS,  90,0.20) 

CALL  PLOT (0, 0.5.0, '3) 

C 

C SET  UP  VERTICAL  AXIS  OF  LASER  PLOT 
C 

1 NARK  • 4 
XN  ■ INAR4/2  > 

RABOVE  ■ .ABOVE ( 1 ) 

ABOVE  ■ R ABOVE *0 . 30*8 
RBElOH  ■ .RELOH(l) 

BELOW  ■ R 3£LOW*0 t 30*8 
AbBh LOW  « ABSF(BELOW) 

SCALAR  > MAXlFlABOVfc.ABHELOW) 

XmIU  ■ 2.0 

CALL  PLOT(-O.l,0.0.3) 

Call  PLOTtn. 0,0. 0,21 
00  10  I ■ 1.4 
YY  * t 

CALL  PLOT  ( O.O.YY.l) 

CALL  PLOT  (.fl.l.YY.l) 

call  plot  ( o.o.yy.i) 
io  Continue 

Call  NUHBjR  (>0.6,  3.95, 

SCALAR  ■ .SCALAR 

call  symbol  (-0.6,  1.95, 

CALL  SYMBOL  <.0.6,  1.58. 

CALL  NUMBER  (-0.6..U.05. 

SCALAR  e -SCALAR/2. 

c 

C PLOT  lasea  record 
C 

XX  s DELT»(NUMPT  ♦ NST) 

LL  = 3 

no  20  1 = NST.NF1NAL 

YAWP  > Y!  I >*0 . 3048/SCALAR.XMllJ 
IF (YAMP  . 3T . 4.0)  12.14 

12  YAMP  ■ 4.) 

00  TO  17 

14  If  (YAMP  ..<T.  0.0)  16,17 

1ft  YAMP  a 0, 

17  CALL  PLOT  (XX, YAMP, LD 

If  <LL  • 2)  175. 170, 171) 

1/0  LL  = LL  - 1 

1/5  If ( PlTONL  Y ,EQ.  0)  20,18 

18  IF (M00( I , N°TS)  .EQ,  0)  19.20 

19  CALL  PLOT  (XX, 4, 5, 3) 

CALL  PLOT  (XX, 4. 0.2) 

CALL  PLOT  (XX. YAMP, 3) 

CALL  PLOT  (XX. VAMP, 2) 

20  XX  s XX  * HELT 
C 

C SET  JP  HORIZONTAL  AXIS 
C INARM  SHOJLD  BE  AN  EVEN  INTEGER 

C 

DELTAX  ■ XXF'MAX/IMAHK 

JHIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
gHKHi  OOPY  FURNISHED  TO  DDC 


0.105,  SCALAR. U. o. 4Ef 6.2) 

0.105,  4H  0,0*0. 4) 
0.21u.5HLASER,90,0,5) 
0.105,  SCALAR. 0. o. 4ff 6.2) 
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CAUL  PLOT  (XXFMAX.-U.5,.)) 

JHANX  a MARK  ♦ 1 

DO  SO  1 al , JM  ARK 

XA  a XXfMAX  - < J - 1 > »DfcL TAX 

Xb  » xa*p.itspace/delt 
Call  plot  <xa,-o.5.2> 

Call  plot  ua,-o.55,1> 
call  plot  <xa.-o.45,1> 

Call  NUMBER  (XA-0. 255.-1), 40.0. 11)5, XB.O.U.AHM.O) 

IM I ,E0. I 4ARK/?  ♦ 1>  25.27 

25  LAll  SYMJOL  <XA  - U. 88, -0,85.0. 14. 15HTJME  I*  SECONDS, 0 . 0 , 15) 
27  CAuL  PLOT  (XA. -0.5.0) 
jO  CONI INUE 

IKnJHHIST.EO.2)  32,55 

set  up  vertical  axis  of  acceleration  plot 


02  Call  plot  <o.,-5.o#»3> 
call  plot  <-0.1.0.0|3) 

Call  PLOT  < 0 • 0 . 0 ■ 0 i 2 ) 

DO  40  la  1.4 
YY  a I 

call  plot  c o.o.yy.d 

CALL  PLOT  (-O.l.YV.D 

Call  plot  ( o.o.yy.d 

40  CONTINUE 

RaBOVE  a .ABOVE ( 2 ) 

ABOVE  a RABOVE*0.3048 
RbELOM  a .BELOW( 2 ) 

9EL0N  a R}ElOW*0.304B 
ABBELOW  a ABSF(BELOH) 

SCALAR  a 4AXlF(AB0VE. ABBELOW) 

Call  NUMBER  (-0,6.  O.V5.  0.105,  SCALAR. 0.0, 4D6. 2) 

scalar  a -scalar 

call  symbol  <-o.6,  1.95.  0.105.  4H  o.u.0.4) 
call  symbol  (-0.6.  i.ss,  o.21u.5maccel»9o.o,5; 

CALL  NUMB-R  ( -0 , 6 , -U , 05 . 0,105.  SCALAR. 0 . 0 . 4FF 6 . 2) 
SCALAR  a -SCALAR/2 • 

XX  = DELTDNUMpT  ♦ NST) 

LL  « 3 

PLOT  acce.eration  RECORD 

00  50  I « NST.NFINAL 
VAMP  a /< I )*0. 3048/SCALAR  « XNIU 
IF (YAMP  . JT.  4.5)  42,44 

42  YAMP  a 4.5 
QO  TO  48 

44  IF  ( YAMP  ..IT,  -0,5)  46.48 

46  YAMP  a .0,5 

48  call  plot  cxx.yamp.lu 

IF (LL-2)  50,49,49 

49  LL  * LL  - 1 

50  XX  * XX  ♦ DELT 

CALL  PLOT<0.0,5.0,-4) 


POSITION  »EN  FOR  NEXT  PLOT 

55  IF (SWHONLr.EQ.l.O)  56.58 

56  CALL  PLOT(n,0,4.5,-3) 
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GO  13  66 

58  If (PLTONLY.EO.O)  60.62 
60  CALl  PLOT(0.0.-l«.0»-4> 

Go  '3  64 

62  IMM.Efl.  l.AND.NT.EU.NTIMES)  64.96 
64  XMXXFMAX  * 4,0 

CALL  PUflf (XF, -0.0,0) 

66  RETuRN 
END 

SUBROUTINE  PLTSPEU  i yp,xp,lphi,hbaR,nop.np. ymax, nomaVplt, 
*XMAXVL, ISEC, I8EL0N. » ABOVE, MS AT, N2. THETA, Nl.NCS.PUf SPACE > 
dimension  yp<4097.2>.xp(4ov7,2),lphi<iu) 

Di HrNSION  H8AH(2> , IbTARl (2J 
type  INTE3ER  THETA 

LPHl  CONTAINS  PLOTTING  LABEL  TO  BE  PLOTTED  ON  TOP  3r  PLOT 
MAXIMUM  on  72  CHARACTERS  PLOTTED  OUT 

TP(l.N)  * ARRAY  OF  V VALUES  OF  THE  N-TH  VARIABLE  T 3 BE  PLOTTED 
XP(I.N)  * ARRAY  OF  X VALUES  OF  THE  N-TH  VARIABLE  TJ  BE  PLOTTED 
NOP  - TJTAL  NUMBER  OF  POINTS 

NP  . STARTING  POINT  IN  ARRAY  FROM  WHICH  TC  START  PLOTTING 
FpAx  « PHYSICAL  LENGTH  OF  ABSlSSA 
XLE>T  ■ L3G  OF  MINIMUM  VALUE  ON  ABSlSSA 
XhAXVL  • .OG  OF  MAXIMUM  VALUE  ON  ABSlSSA 
2U0  FORMATI/. 2X..ISTARTI1)  ** . 1 3, 5X, *NOPP  **. 1 9, 5X, * I STaRTC 2)  «*,I3. 

2 5X.*N0P  »«  » I 9 » / ) 

5V0  FORmAT(3X,«caSTplT«  *.I5> 

CONVTM  > 3,3048 
CONVETM  ■ CONVTM**2 
CONV  * l . /LOGF (10  , ) 

AlHtTA  « THETA 
If (THETA. 3T.180)  525 ■ 590 
529  THETA  > 3JT  - THETA 

550  Do  3,  N«l,2 

2 ISTaRT(N>  = 1 

3 CONTINUE 

CA».L  NEwPiN(t) 

6 FMAx  * 

CONVERT  Data  TO  METRIC  SYSTEM 

HBARU)  « CONVTMPHBARU) 

HbAR(2)  x CONVTM* HBAR(2) 

DO  10  I * l.NOP 
VP(l.l)  * CONVETM*YP( 1 ,1) 

10  YPII.2)  » CONVETM*VP(  1,2) 

YMAx  ■ CONVETH*VMAX 

DETERMINE  INITIAL  V AND  Z VALUES  FOR  PLOTS,  CALCULATE  LUG  VALUES 


XMAXVi.  • },9 
XLE-T  ■ -2.9 
XLS^TIO  ■ 10.**XLEFT 
SCAut  ■ 5, 

DU  103  N ■ 1,2 
DO  191  1*1, NOP 

If<X»<t,N>  ,lT.  XLEfTIO)  93, 6U 
50  ISTaRT(N>«  1*1 
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IBIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
ISO#  COPY  PUKKISHED  TO  DDC 


00  TO  too 
6tl  XPl  a XP(I.N) 

XPL  ■ UOG?<  XPU 
XPl  1 ,N>»CJNV*XPL 

too  continue 

103  CO NT  I N J£ 

IMn2.Eq.I>  104.10« 

104  IMhSXT.GT.  ISTARTIin  105.11)6 
lot*  ISTaRTII)  » M$AT 

106  ISURTl  * I START ( 1 ) 

SfeT  JP  0OTH  AXES,  PLOT  AND  LABEL  VERTICAL  AXIS 

DO  108  I « ISTART1.N0P 
U (XP(I.l)  .QT.  XMAXVL)  107,108 

107  NQPP  a 1-1 
GO  TO  10V 

108  CONTINUE 
NOPP  « NO >i 

XMAXLVL  » XHAXVL-XLfcFT 
XMlN  « F«AX/XMAXLVL 
GO  TO  110 

1U9  XNAXLVL  * XP(NOPP.l)  * XLEFT 
XHIN  * TMAX/XNAXLVL 
HO  YHAX  ■ LOIF(YNAX) 

YMAX  a COxVaYMAX 
YNAx  a YMAX«1, 

JYMAX  a Y AAX 
YT  OP  a JYAAX 
YBOTTOM  » YTOP-SCAlt 
YOOTTOHR  ■ 1 0 , »a YBOTTOM 
CALL  PLOT  ( 0.0,  -.0,1.3) 

Caul  plot  ( 0.0,  1,0.2) 

DO  111  I a 1,9 
Y « I 

Call  plot  <*o.l,  Y.l) 

CALL  PLOT  C 0.0,  Y.l) 

Call  plot  < o.o.y*i,,i> 
m continue 

CALL  PLOT  (-0,1  . Y*l.  , 1 ) 

call  SYM8JL  <0. 3. 10,1,0. 14. LPHl.O. 0,72) 

Call  symbol  (9.08,9,57, 0,14, 0.0.0, -1> 

CALL  SYMBOL  (9.15.9,50,0.14,26H  ■ LASER  ( SNH  ■ M,  >,0.0,26) 

CALL  NUMBER  (11,07,  9 . 50 , 0 , 14 , HBAH < 1 ) , O . 0 , 4HF 5 . 2 > 

CALL  SYMBOL  (9.08,9,32,0,14,4,0,0.-1) 

CALL  SYMBX  <9.l5,9,2b,r).14,3UH  • A/C  Utsp,  (STD  « M,  ),0.fl, 

*30) 

CALL  NUMBER  (11,60,  9 . 2b, 0 . 14 , HBAR ( 2 > , 0 . 0 . 4HFS . 2 > 

CALL  SYMBDL(9,13.9.a0.0.14.8HTMfcTA  a ,0.0,0) 

CALL  NUMBER(10.15.9,00,U,14,ATHETA,0.0,4Hr5.1) 

YY  a 9,9b 
FlABEL  a /TOP 
DO  116  I a 1,6 

CALL  NUMBER  (-O.8,YY,O,lO5.FLA0fcL.O.O,4HF6,2) 

FLApEL  a CLA8EL  - SCALE/5. 

YY  = YY  - 2,0 

IMI  .EQ.  4)  119.116 

lib  CALL  SYMBX  (-0.8, 3, 6,0.173, 22HL0G ( SPECTRA ) M**2-5tC, 9u , 0 , 22 ) 

116  CON) INUE 
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TRANSFORM  DEPENDENT  VARIABLE  TO  log  BASE  10,  CONVERT  COORDINATES 
lo  INCHES. 

JJ  = NOPP 

DO  121  N • 1.2 

If (N.E0.2)  117.119 

n "/  jj  = nop 

ll»  ISTaRTI  « ISTARTIN) 

00  120  I • [START!. JJ 
YPL  • VP ( I » N ) 

If IyPL.LE.TBOTTOMR)  1002.1000 
1000  YPL  » L0G-*<YPL> 

YPd.N)  « C0NV* YPL  « YBOTTOM 
YMI.N)  = YP(  I ,N)/SCALE*10. 

GO  TO  100* 

1002  YHi.N)  * 0, 

1 0 0 A XP(I.N)  * XP(I,N)  * XLEFT 

120  Xpd.N)  * XPfI,N»«XMIN 

121  CONTINUE 

FMll'  a FMAX/2, 

J1  « 0 

Nfc  NL>  ■ NO»P  - NP  ♦ 1 
Jj  = NEND 

If(  THETA. LE. 90. ANU. 1S1ART(1) .EO.l)  122.210 

122  ISTaRT<1>  « ?. 

210  PRINT  200,  ISTaRTUj.HORP,  I START  { 2 ).  NOP 

PLOl  SPECTRA  (LOG-LOG) 

INZ=N0P/2S6 

ANZ=INZ 

00  140  N » 1.2 
L e -1 

If (n.EQ.2)  123.125 

12J  J 1 -4 
Jj  : NOP 

125  IST4RT1  = ISTART(N) 

DO  135  I ( ISTARTl  .JJ 
NN  r 1 

If (N.EQ.2)  127.129 
127  NN  s JJ  - I * ISTART( 2) 

129  IMYPINN.V)  .LE.  O.f  131,133 
131  L * 0 

GO  TO  135 

133  If (N.EQ.2)  134,500 

LIMIT  THE  NUMBER  OF  HIGH  FREO  PTS  PLOTTED  (LASER) 

5(|0  PRNUH1«NN 

PRNuH2«PRVUrtl/ANZ 
I pRl bPRNO  42 
PRNuH3» 1 P R1 

PRNuH4»pRvUH2-PRNUHO 
If (NN.LE.125.0R.PRNUM4.EQ.O,)  134,135 
lj4  CALL  SYMBOL  ( XP ( NN, N ) , YP (NN, N ) , U , 07 , J I , U . 0 . L > 

LASIPI.T  * NN 
135  L • L - 1 
140  CONTINUE 

PLOT  AND  .ABEL  HORIZONTAL  AXIS,  POSITION  PEN  f OR  NEXT  PLOT 
If  XLEFT  JR  XMAXVL  IS  CHANGED,  LIMITS  FOR  DO  LOOPS  MAY  CHANGE 


BIS  PAGE  IS  BEST 

YBPii  COPY  FURNISHED  10  DDQ 
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XINC  » FKAX*0 ,5/XMAXLVL 
call  PLon-o.i,  o,  i 3 ) 

call  plot  <xinc.o..*> 

x « o. 

DU  141  I • 1.5 
X ■ X ♦ XINC 

cali  plot  < x . -o.i  # l > 

Call  plot  ( x . o,  * l > 

CALL  Plot  (X  ♦ XINC,  U..1) 

141  CONTINUE 

call  plot  <x  ♦ xinc.  -o.i.  i> 

X » -0.3 

XX  = 0, 

DO  150  I • 1,7 
FREQ  * XX/XMIN*XLEFI 

Call  NUMBER  <X, -0.20.0. 105. FREO.0.O.4HF6. 3) 

IP  < XX  ,GT.  FMIU-1.5  .AND.  XX  .IE.  FMIU)  145.14V 
145  CALL  SYMBOL  (FMlD-1, 47. -0.5,0. 175, 20HLOli<FRECUENCY>  M2  ,0.0,201 
14V  XX  s XX  * XINC 
150  X » X * XINC 

IMnOWAVP.T.EQ.I)  155,157 
155  CALL  PLOTIFMAX  ♦ 4.,0.0.-3> 

Go  '0  180 

IS'/  If  ( I SEC  , FQ . 1 > 160.170 

160  IMM.EO.l)  162,165 

162  CALL  Plot  (FMAX  ♦ 5,0,0.0,-3> 

GO  TO  180 

165  XflN  » N05»PLTSPACE  ♦ 5.0 
CALL  PLOT  (XF IN, 0.0, -3) 

GO  TO  180 

I/O  NPOINTS  = l IAB0VE*I«EL0W)#4 
INCHES  ■ VPOINTS/4 
1 1 ST  AN  ■ I NCHES*9 
DJSIAN  ■ I I STAN 

call  plot  <fmax*distan,o.o.-3> 

1BO  RETURN 
END 

SUBROUTINE'  PLOTDFIDM, 1 BELOW, I ABOVE . M I N. MAX. LPHl. HPlOTuF ) 

RETURN 

END 

SUBROUTINE  ROTaPE 

CONMON/RUTAPEB/LFILfcl  7I.XFILE.  IHLE.MFILE 
IMmFILE-<FILE.EO.O>  30,10 
10  K1»MFILE-<FILE 
00  70  1*1, <1 
20  CALL  SK I Prll  LE  (5  > 

Kf I L E*MF I .E  * 1 
30  CALL  UNpC  < (0,0,01 
RETURN 
END 

SUBROUT  I NEi  UNPCKUSWP,  IKEC.NFLAU) 

DIMENSION  INSK8.B),  INS ( 264  ) , VSCOPE { 9 ) , 1 NS2 ( 8 ) ,NRA  ( llOOu  ) . 
lNRA[i(  167 ) 

COMnON/uN>CKB/RAD(160),CHAN(7),AlNSS(8), I NCOTE ( 8 ) » NF 1 LE , 
lNUAy, MONTH, NYR.NMR.MIN.ALT.UDOT.NH SCOPE. NVSCCPE.VSjALE.jS WEEP 
COMMON/B/ J< 550 0.2 1 ,»  NUHTF(40»7,2) 

EOUlVALENOE  <NRA.U»)NUHTF> 

data  (M1«17B), <M2*7/77R), (M3*1B), (M4-3T00B), (M5«lO00B) 

DaTa  ( M6* 77777777777760  0 OB  1, (MSK1*777740 000  OB), 

C (MSk 2*777 7740 OOOB), I MSK3*377B), (MSK4*14000000000B>, 
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DIMS*  5 = 77/  7 7700  00000U0PB) 

Da  r M (He*  I'H) . < ( VSCUPtl  I>.l«l,9)*l.,2t.5.,10..20..50.»lu0,,200,, 
C5u0  . ) 

U<  iREC.Ej. 0)2000. 1020 
2000  CONTINUE 

B U F . EH  IN(S, L) (n»A(1) .NRAI2)  ) 

110  IF (UNIT, 5>Jin, 30.30*40 
*0  PRUT  50 

SO  FORMAT!  • PARITY  ERROR  IN  FIRST  RECORD*) 

JU  DO  l J»l.J 
DO  1 1*1.9 

MTS=NRA(  J)/2**(  ( 1-1  >*6) 

INSl ( I , J) »RTS.ARD,M1 

1 CONTINUE 

NHlE«INSII4.1)  %NDAY*INSl(H.l)»tO*INSl<7,l) 
MONTH*1nS1<6,1>*10*)NS1I5,1)  SNTR*7 0*1 NS1( 3 , 1 > 

NRRs INSl (2.1) *10* INSl 1 1.1) 

MlNslNSU 9.2>*10*IRS1(7,2> 

ALT=INSl<  5,2>*tO*INSl(5,2) 

DDOI»INS1(4.?)«10.*1NSK3.2) 

NHSi,OPE»l  YSl(2.2) 

NVSC0PE*I ySI (1,2) 

vscale*vs:ope(nvscore> 

PRINT  100.NF ILK, NDAY, MONTH, NYK 

100  FORMAT (lX.jHF  Ii.E  , I2.2X, 12.1H/, J2.1H/, 12.//) 

PRINT  101.NHR.MIN 

101  FORmATIIX, 7HT1ME*  « 1 2 , 3HHR  .2X.12.4HMIN  ,12.//) 

PRINT  102.DDOT , NV SCOPE , NH SCOPE .ALT 

102  FORMAT(2X,10HDFL*T  001«,F4.1,2X,11HVErT  SCALE*. 

C1X.  12.2X.  IOHHOR  SCALE*. IX,  1 2 , 2X,  4H AU T» . F5 . 1 , / / ) 

IRCsO 

return 

1020  CONTINUE^  IF(IRC.EQ.IREC)  75,52 

IRC:  I REt^  52  IRC2  = IREC  - IRC 

DO  n IBUrPl, IRC2 

55  BUFFER  IN(5,1) (NRA(l).NRA(llflOO)  > 

60  IF (uNIT, 5)60,7 0,65, 60 
BO  PRINT  90 

YO  FORMAT ( • PARITY  ERROR  IN  SHEEP  DATA*) 

GO  TO  7o 
65  CONTINUE 
/0  CONTINUE 
H CONTINUE 
/5  CONTINUE 

MM* ( IREC-l)*ISWEEP*JSHP 
Jj*M0D(MM,/56) 

DO  301  J*  t , NFLAG 
NUM=< ISwP'l ) *167* J 
1h0R0*(NU4*3)/4 
HSM 1 F T * i 4 * 1 WORD -NUM ) *12 
NRAUI J)»NTA( IWOR0)/2**NSHIFT 
M)S:NRAD<  J> 

NRAi.(  J)*HTS,  AND. M2 
3U1  CONTINUE 

IF IJj.EQ. 1)302,109 
302  1*0 

DU  3 J*l. 132 
1*1*1 

HTS  = NRAD<  J)/?**ll 
INSl 1 )«MTJ.ANU.M3 
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10 


210 

11 


1 A A 

14 

15 

16 

17 

1H 

IV 

20 


22 

24 

26 

20 

2 V 
01 
02 

03 

1 03 
109 


12 

13 

114 


US 


1 ■ I *1 

Ml  Sc  NR AD  I J 1 / 2 • • 5 
I NS ( 1 > »mT  5 , AND . M3 

continue 
1 1*7 

no  u k*i,b 
1 NS2( K 1*0 
DU  ’0  1*1,8 
X 1 ■ 1 1*1 

1NS2(K)«I NS(  111*2**18-1 1*1  NS? <K) 
CuNT INUE 
DO  210  l = >,32 
U«ll*l 

INS'2(K)*1  4S(  tl)*2*«t  I*l)*lNS2lK> 

Continue 
coni INUE 
DO  20  K*1 , 8 

INCuDE(K>»MTS.AND,MSK3 
MTScINS2( <) 

IX U.tQ.2.0R.K.EQ.5)133,14 
INS2(K)*MTS,ANI),MSK2 
GO  TO  15 

INS2(K)«MTS. AND , m ski 
NSU.NaMTS.AND.MSK4 
If  USIGN.NF.  0)16,17 
Ml  Sc  I NS2  < 0-1 
INS2(K)*MTS.0R.MSK5 
CONTINUE 

If (2. fcQ. 2, OR. K. £0,5)18,19 
AlNSS<Kl*|NS?<K)/2**29*.127ft,7 
GO  I 0 2 11 

AlNSS(K)  = I'tS?<K>/2«»2V«180. 

CONI INUE 

AINsS<5)*AlNSS<5>*2. 

If  < a INSS<  3) ,LT.O . >22,24 
AlNSS<3)=AINSS(3)*36fl. 

If  < A INSS<  *) .LT.O. 126,28 
AlNbS<4)aAlNSS(*)*360. 

If  (A  INSS<  11  .l_T  . 0.  129,31 
A|NSS(1)=AINSS(1)*360. 

If  < a INSS( S) .LT.n. >32,33 
*INSS<8>=AiNSS(8)*3O0. 

CONI INUE 

PRINT  103, ( A I NSS ( I ) I 1=1.8) 

FjRmA  T <2X, 4HINS*. 8(2X.Fa,l) ) 

Du  114  11*1,167 

MIScNRADT 111 

Ml  Si *MTS. AND.M4 

MlS^iMTS.AND.MC 

NHAi  (U)*MTS1/2»MTS2 

Ml ScNRAD ( 1 1 1 

NSU.NaMTS.AN0.M5 

If  INSIGN, NP. (1)12,13 

MlSl*MTS-l 

NR  A i ( Il>aiTS’  .OR. ,4ft 

CONI INUE 


H INUE 
115  J6  a 1 , 7 

N( J6)=4RA0< J6)«l. 953125 

THIS  PA®  IS  BBS!  OVIAilTl  l MCTlCMa* 

jgQH  COPY  FUKfrlSHED  TO  DD.C.  ^ 
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DU  107  I«t.l<S0 

RADI  I )»NRA(><  l*7l«VSU*l.e/51.Z 
1 U 7 continue 
RETuRN 


END 

SUdHOUT I N • PHASE<NTOl ,UP) 


COMmON/A/Y.Z 

U 1 Me  MS  I ON  Y(8192>#Z<0192) , II YPE<2) 

Data  < !TY»f«<5HlASER.5HACCEl> 

9«u  format  (»3  atest  ■•ra.i,*  rtest  •»f6.i<*  STEST  **F'/,l,//> 
9u 1 Format  <«J  NUMHtR  Of  PHASE  SHIFTS  FUR  aA5«a  data  « •»I4./» 

• 6X,aM4X  >HASh  SHIFT  = • , F6 , ? » 5X , «M I N PHASE  SHIFT  « *,Fo,2,/, 

• 6X, AUIff jRENCFS  GREATER  THAN  RTEST  « * » 1 4/ > 

20,10. JO 


If  (UP 

• 

13.  ) 

10 

Al  EST 

• 

4.0 

RlfeST 

* 

2. A 

S1EST 

a 

O’ 

U1  EST 

a 

32.2 

TIEbT 

9 

64.4 

B1EST 

a 

32.2 

GO  10 

«0 

20 

AlfcbT 

a 

3.5 

R1ES1 

a 

2.0 

S1EST 

a 

Q> 

U1EST 

a 

32.2 

TlfcsT 

a 

04.4 

BleST 

a 

32.2 

Gli  10 

40 

60 

A1EST 

a 

5.0 

R1EST 

a 

3.? 

SlEsT 

a 

06 

U1EST 

8 

32.2 

T1  fcST 

a 

04.4 

BlEfT 

8 

32.2 

All 

PRINT 

90  0, 

ATEST.RCfcST, STEST 


Ml  0 1 = NTJT  - 1 

NL  ASPS  ■ YRAllPS  * NLASDlF  * NRAUD l F * U 
RYMaX  * RZMAX  = 0. 

Rymiv  • RZM I N = 1000. 

60  DOT  * DCZ  = 0. 

DO  ISO  I * I.mTOT 

V « Y<  I ) 

H * Z ( I ) 

Y ( I ) s Y C I > ♦ DC Y 


U L * ACCELERATION,  NO  PHASE  SHIFT  IS  NEEDED  ANE  DCZ  « 0. 
Di;Z  s il , 


zii)  * zl l ) ♦ nez 

AY  - Y( I) 

AZ  = Z(I> 

dY  s Y C I ♦ 1 ) * I1CY 
BZ  * Z» 1*1>  * DCZ 
Ry  = A3Sf (AY  - BY) 

RZ  * A9Sf (*Z  - BZ) 

SY  s AY  - HY 

SZ  A AZ  - «z 

If ( I COUNT  .EO.  2) 
/»  If (MY  .Li . ATEST) 


1X0,70 
MU,  75 


THIS  PAGE 
FROM  COP* 


IS  BEST  QUALITY  FRkCIICAW® 
jUESNISHEDrQSSfi ^ 
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oncinooocicio 


/•»  DUY  ■ OCY  * (JP*  ( RY/oiY ) 
NLAbPS  a MU  ASPS  ♦ 1 


u <ry 

,GT  , 

RYMAX ) 

76,  >7 

/ 6 

RYM<,X 

a RY 

GO  10 

80 

/ 7 

If  I ft  Y 

.LT, 

RYMIN) 

73, H9 

/b 

RYM  i N 

i R Y 

60 

If  (HZ 

.lt. 

BTEST) 

loO . 90 

VO 

DUZ  a 

DCZ 

♦ RZ* (RZ/SZ  > 

NHAi.PS 

; a mPADPS  ♦ 

1 

If  I ft z 

• GT , 

rzmax) 

95,100 

VS 

R/MAX 

a RZ 

Go  TO 

160 

luo 

If  (HZ 

.LT. 

RZMIN) 

105.160 

105 

rzmin 

a H2 

Gil  1 0 

160 

120 

If  ( ft Y 

. G T , 

RTEST 

.AND,  RY 

.LT. 

STEST) 

130  i 

-140 

160 

DU  Y a 

DCY 

♦ ry*ry/sy 

NlASDI 

K a 

NLASDIF 

♦ 1 

1*0 

If  l»Z 

. G T , 

UTEST 

• A N i J # HI 

.LT. 

TTEST) 

150, 

>160 

laO  Du  Z a DCZ  ♦ RZ*RZ/SZ 
NPAuDir  * NRADDIK  ♦ 1 

160  continue 

Y(NTOT)  * Y ( N TOT  ) ♦ DUY 
C ZINTOT)  * ZINTOT)  * ouz 

IMICOUnT  .DO.  2 ) 200.1  86 

160  IUDONT  ■ ICOUNT  * 1 
Go  *0  69 

200  PR  I NT  991.  ITYPE»1).NLASPS.NYNAX,HYMIN.NL*SDU 
PRINT  991.  ITYPE<2),NHADPS,RZMAX.RZMIN»NRaDDIK 

re  Turn 

END 

SUBROUTINE  FOURIER  (A.S.H.tfS) 

THIS  ROUTINE  PERFORMS  AN  ANALYSIS  OK  2*»«  POINTS  NY  K IHST  DOING 
AN  ANALYSIS  ok  2»*N/2  CONPLfcX  POINTS  AND  THEN  ARRAMOiNG  THE  RESULTS 


ARGUMENTS 

1.  A - REAL  DATA  ARRAY  - OK  DIMENSION  2**m  « 2 
27  S - SI  M/COS  TARuc  * DIMENSION  2**<M-3> 

3.  M - EXPONENT  OK  2 - SIZE  OK  HEAL  ARRAY 

4.  IKS  - -1  KOR  FIRST  TIME,  -2  YMENAFTER 

DIMENSION  * < 1 ) . S( 1 ) 

N ■ 2**(M.1 ) 

CALL  HARMJNIA.S.M-l!  IKS,  lf'EHR  ) 

C MERGE  2 N. POINT  ANALYSIS  INTO  1 2N-P0INI  ANALYSIS 

NMAlK  a N /? 

NT  NO  a N * 2 * 4 

X * XO  a :0SF<3.14159'26536/FLUATF<N)) 

Y • YO  ■ SINK (3, 1419926536 /KLUATFIN)  > 

00  1009  K2  a 4.N.2 

K1  = K2  - 1 

N2  * NTUO  • K2 

N1  a N2  - 1 

BKl  a A ( K 1 > ♦ A ( N1 ) 

BK2  a A < K 2 ) * A(N2) 

BN1  a A ( X 2 ) * A ( N2 ) 

BN2  a A(X1>  • A(N1) 

XBNl  a X*j|Nl 
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THIS  PAGE 

KBUM  COF*  yUBJUSHBU  w 


XMi'  « X»JN2 
yfa*«i  * y*j*i 
YbN2  » Y*JN2 

AK.I  * .5  • ( HK t ♦ XBN1  - YHN2) 
AlKi)  ■ . S * ( *BK2  * XBN2  * Y0N1J 


Al2) >*.S 


A Hi)  ■ .i  * ( -BK2  * XBN2  * yBND 
A(Nj)  « . J • (HK1  - X8N1  ♦ YBN2) 

AIN*>  ■ .}  *(BK2  * XBN2  * VBNl) 

0 » X*X0  • Y*YD 
Y ■ Y»X0  * X*Y0 
X ■ 0 

C04PLEX  E.FM6NT  A(N> 

A(2*N*1>  i ( A ( 1 ) - Al2>)*,3 
A A 2*  >1*2 ) * 0.0 
COMPlEX  E.C'MENT  A<0J 
All)  « ,!»IAI1)«A(?I) 

A ( 2 > « 0. 3 

COMPLEX  E.EMENT  A(N/2) 

A « N* 1 ) ■ *(N*1> 

A I N*2 ) « 4(N*2> 

RETuRN 

End 

SUBkOUTIN;.  HARM0N(A,S.M. I FS, IF  ERR) 

DINi-NSION  A(l).stt) 

harm,  dne -dimensional  basic  puhtran  version,  j.w.cuoley  harm  ooi 
MODIFIED  ro  RUN  ON  CDC  -SlUD, 

haRH  002 
HARM  006 

DUES  EITHER  FOURIER  SYNTHESIS. I .E. .COMPUTES  COMPLEX  FoUKlER  SERIESHaHM  009 

OIVcN  A VECTOR  OF  N COMPLEX  FOURIER  AMPL I TUDtS, OR.  UIVEn  A VECTOR  HARM  010 

of  complex  data  x Dues  fourieh  analysis,  computing  amplitudes.  harm  oil 

A IS  A COMPLEX  VECTOR  01  LENGTH  N«2**M  COMPLEX  NOS,  OR  **N  REAL  HARM  012 

NUHhERS,  A IS  TO  BE  SET  BY  USER.  HARM  013 

M IS  AN  INTEGER  0.LT.M.LE.13,  SET  BY  USER.  HARM  014 

S IS  A VECTOR  S(J)»  S l N ( 2*PI  • J/NP  ),  J«1 . 2.  . . . . , NPM«1,  HARM  015 

COMPUTED  9Y  PROGRAM.  HARM  616 

IPS  IS  A PARAMETER  TO  BE  SET  BY  USER  AS  POLLCHS-  HARM  017 

USsO  TO  SET  NP«2**M  AND  SET  UP  SINE  TABLE  S.  HARM  018 

IFSil  TO  SET  NsNP»2**M.  SET  UP  SIN  TABLE.  ANC  DO  FJURlEH  HARM  019 

SYNTHESIS,  REPLACING  THE  VECTOR  A BY  HARM  020 

HARM  021 

XIJ>«  SUM  OVER  X*0»N-1  UP  A|K>*EXP(2*Pl*l/N)**tJ*K)t  HARM  022 

J»0,M-t,  THERE  I»SQNT(-1>  HARM  023 

THE  X«S  ARE  STORED  WITH  PE  X(J)  IN  CELL  2*J*1  HARM  024 

AnO  IH  X(j)  IN  CELL  I.'* J*2  FUR  J«0,1.2«..'.N"1,  HARM  025 

THE  AXS  ARE  STORED  IN  THE  SAME  MANNER.  harm  026 


IF  3s-L  TD  SET  N»NP*2**M.SET  UP  SIN  TABLE,  AND  00  FuUhIER 

analysis,  taking  the  input  vectur  a as  x and 

REPLACING  it  by  THE  A SATISFYING  THE  ABOVE  FOURIER  SERiES, 

IF  Ss *2  TO  DO  FOURIEH  SYNTHESIS  UNLY,  WITH  A FRE-COMPUTEU  S. 
IfSs-2  TO  no  FOURIEH  ANALYSIS  ONLY,  WITH  A PRE-cOM»UTED  S, 
IFEhR  IS  SET  BY  PHOGRAM  TO- 
•0  IF  NO  iRROR  DETECTED. 

■1  IF  M IS  OUT  OF  RANGE.,  OR.  HHEN  lFS**2,-2.  THE 
PRE-COMPUTED  S TABLE  IS  NOT  LARUE  ENOUGH. 

>•1  WHEN  I FS  »*1,-1»  MEANS  ONE  IS  RECOMPUTING  S TaSL^ 

unnecessarily. 

NOTE-  as  stated  ABOVE.  THE  MAXIMUM  VALUE  OF  P FOR  THIS  PROGRAM 
ON  THE  IB*  7094  IS  13-  ON  360  MACHINES  HAVING  QHEaTER  S, OH-AGE 


HARM  027 
HARM  028 
HARM  029 
HARM  030 
HARM  031 
HARM  032 
HARM  033 
HARM  034 
HARM  035 
HARM  036 
HARM  037 
HARM  038 
HARM  039 
HARM  040 
HARM  041 
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no  o no  onnnonnn 


IMIS  pass  IS  BEST 

Wl! 


CAPACITY*  ONP  SHOULD  CHANGE  THIS  LIMIT  BY  MeFL  AC  I N 3 13  < N 
STATEMENT  3 BELOW  BY  L0U2  N,  WHERE  N IS  THE  FA /.  NO,  OF 

COMPLEX  NJMBERS  one  can  STOHE  in  high-speed  come, 

IF  HE  CAPACITY  OF  HANM  ts  To  BE  INCREASED,  ONE  MUST 
ALSU  ADD  40RE  DO  STATEMENTS  TO  THE  BINARY  SORT  KDUTiNfc 

following  statement  24  anu  change  the  equivalence  statements 
for  the  <*$. 

dimension  k < 12 » 

EQUIVALENT  (K< U) ,R1) , (K(1U) ,K2) . (K(9) ,K3), (KC8), <41 ,<M />,K5) 
EQUIVALENT  (K(6).KO). (K(5).K/>, (KC4>,Kd), (K<3),<V>, (K(2),K10> 
EUUlVALENOt  (K(1),K11),(K(1),N2) 

If  IP  >2,2, 5 
3 If(M-il)  S.5.2 
2 If  EHR*1 

1 return 
9 If  ERR*0 
N«2**M 

I f < X ABSF  < I F S > - 1 > 2U0.200.1U 

WE  ARE  DOING  TRANSFuRH  ONLY,  SEE  |F  PPE-C()PPL  TED 
S Table  is  SUFFICIENTLY  LARUE 
1U  If { N-NP  >?U.20,12 
12  If  ERR»t 
GO  TO  200 

SCRAMBLE  A,  BY  SANDEXS  NETHOO 
2U  KU)»2»N 
DO  22  L»2,M 
22  K<L)«K(L-t>/2 
DO  24  L*M« 10 
24  K(L+1>«2 

NOTE  EQUIVALENCE  OF  KL  AND  K(l4*l> 
binary  sort- 
I j«2 
Jl»2 

2 DO  30  J2*  JT  ,K2,K1 
DO  .<0  J3*J2,K3,K2 
DO  >0  J4*  J3,  K4 , K3 
DO  .'0  J5*J4,KS,K4 
DO  .'O  J6*JV,K6,K5 
Du  sfl  JfaJ6,K/.K6 
Du  30  Jrt3 J7.KO.K7 
DO  30  J9»J4.K4,K8 
DO  30  JlOUR.Kl0.K9 
DU  30  JI* Ji 0.K11.K10 
If ( I J-JI >2H, 30,30 
20  T»  A ( | J-i  ) 

A( I j-l)«A( J|-l> 

A( Ji-1)«T 
T»A ( I J ) 

A{ iJ)*A( JI > 

*(J| >«T 

30  I J» | J*2 
Jl» jl*2 

If  <*l-Jl>Ji ,29,25 

31  If ( lr$>J2, >, 36 

C 00 1 jG  FOURIER  ANALTSIS.SO  OIV.  BY  N AND  CONJLUATE. 

32  Fn  s FLOAT, <N) 

Do  S4  J*1,N 

A t 2- I -1 > » A(2*I-l)/FN 
34  A(2*I)>-A(2«T)/FN 


HARM  #42 
HARM  043 
HARM  044 
HARM  045 
HARM  046 
HARM  047 
HARM  046 
HARM  049 


HARM  095 

HARM  057 
HARM  05S 
HARM  099 
HARM  0*0 

HARM  062 
HARM  063 
HARM  064 
HARM  065 
HARM  066 
HARM  067 
HARM  0 AO 
HARM  069 
HARM  070 

HARM  072 
HARM  073 
HARM  074 
HARM  075 


HARM  078 
HARM  079 
HARM  080 
HARM  001 
HARM  002 
HARM  003 
HARM  004 
HARM  005 

HARM  009 
HARM  090 
HARM  0*1 
HARM  0*2 
HARM  0*3 
HARM  0*4 
HARM  0*5 
HARM  096 


HARM  0*7 
HARM  0*8 

HARM  100 
HAHH  101 
HARM  102 
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special  case-  l«i 

HARM 

103 

DO  <tO 

HAMM 

10« 

T * At2«l*1 > 

harm 

105 

*t2»l-t)  *T  ♦ A ( 2* I»i ) 

HARM 

106 

A»2«I*l)*r-A(2*I*l» 

harm 

107 

T»A  t 2« I ) 

harm 

100 

A(2*l)  ■ T ♦ A(2*|«if) 

HARM 

109 

A { 2» I *2 ) = T - A ( 2* I *2 ) 

HA«M 

no 

II (H-l)  2,)  ,50 

HARM 

in 

SET  FOR  L*2 

MA»H 

112 

LEXP1>2 

HARM 

113 

LfcXPl»2*«(L-l  > 

HA«M 

114 

Ltxpca 

harm 

115 

LfeXP»2*«(.*i) 

haRh 

116 

NPLs  2**Hr 

haRm 

117 

NPL  » NP*  2**-L 

HARM 

118 

DO  1 31)  Ls2»H 

SPECIAL  CASE-  Js 0 

HARM 

120 

DO  bO  I b2, M2 • LEXP 

HARM 

121 

U«l  * LEXPl 

HARM 

122 

I if  * 1 1 ♦ LEAPl 

HARM 

123 

Ip  =I2*LtXPl 

HARM 

124 

T«A( 1-1) 

HARM 

125 

aa  l-i > ■ r *a 1 1 2*1 > 

HARM 

126 

A(l,*l)  « T» A ( I 2-1 ) 

HARM 

127 

T * * « 1 > 

HARM 

128 

a < i > « t*a( i?) 

HARM 

129 

A( 12)  * TaA ( 1 2 ) 

harm 

130 

T«  -A(I3> 

HAHM 

131 

T)  * A ( | 3- ) ) 

HARM 

132 

At  1 3*1)  * A(U-i)  - T 

harm 

133 

At  IP  ) • At  11  > - Tl 

HAHM 

134 

Atll-l)  « At  11-1)  *T 

HARM 

135 

At  1 J ) « At  11  > *rt 

HARM 

136 

II  (L-2)  12(1,120,90 

HARM 

137 

KLAST«')2-.EXP 

HARM 

138 

Jj«NPL 

HARM 

139 

DO  110  J»  A . LEXP1 , 2 

HARM 

140 

NPJ.,«NT-JJ 

HARM 

141 

UR«MNPJJ> 

harm 

142 

U1«S«JJ) 

HARM 

143 

ILAST«J*K.AST 

HARM 

144 

DO  1 0 II  Is  J,  ILAST.LbXP 

HARM 

145 

1 1 ■ 1 *LEXP 1 

HARM 

146 

12*I1*LEX»1 

HARM 

147 

IP* 1 2*l£X>3 

harm 

148 

T« A ( I 2-1 ) «UR- A t 1 2 ) *0  J 

harm 

149 

T J »*■  ( I 2-1  )*UI*A<  1 2 > "UK 

harm 

150 

At  1 2-1 ) «A ( I-D-T 

HARM 

151 

At  12  ) s A ( 1 ) - TI 

HARM 

152 

At  1-1)  « A ( 1 -i ) *T 

HARM 

153 

All)  • A t I > ♦ T T 

HARM 

154 

T«-»<  I3-1)*0!-At  13) *UR 

HARM 

1S5 

T 1 »A  < 13-1)»UR-At I3»«UI 

HARM 

156 

At Ip-1)«A( U-l)-T 

B1HM 

157 

A(l.>)  s»<  It  ) -T I 

HARM 

158 

At lj*l)«A( i 1 — 1 ) ♦ T 

HARM 

15V 

A t 1 1 ) —At  1 1 ) *T) 

HARM 

160 

End  dp  i .nop 

harm 

161 

Jj»..J*NPL 



162 
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fflis  PA® 
IBQtf  COPY 


IS  BEST  QUALITY  FRACttCAKUl 
J^PWISifflp  10  XffiQ 


c 

end  of  j .ijop 

HARM 

163 

Iff  u 

LfcXPl«2«tjXPl 

HARR 

164 

LEXP  a ?».1XP 

harm 

165 

t JO 

NHU*v|Pl./2 

HARM 

166 

c 

ENO  OP  L .OOP 

If  ( IF$)145.2,1 

HARM 

167 

cc 

001X3  POlHifcR  ANALYSIS.  REPLACE  A BY  CONJUGATE. 

HARM 

169 

H5 

190 

Dl)  1 50  1*1. N 

A t 2« I ) a- A ( 2» I ) 

GO  13  i 

HARM 

1»0 

c 

RETuRN 

HARM 

173 

c 

maxi-  TABLE  op  S<J>*SlN<2*Pl*J/NH>,J*1.2..,..AT-l,NTaNH/a 

harm 

174 

2UU 

NP*(s 

HARM 

175 

Mpaf 

HARM 

176 

N 1 * X / 4 

harm 

7.77 

Mlaf-2 

harm 

176 

It  <mT ) 26J.260.205 

harm 

179 

2U5 

ThETA*. 7853081634 

harm 

160 

C 

THE! A*PI/2**(L+1>  fOR  Lai 

JSTtP  a NT 

harm 

lSl 

C 

JSTbP  a 2** ( MT-Lfl  > TOR  L*1 

harm 

163 

JO I i * NT/9 

HARM 

164 

C 

JO 1 1 * 2*4(MT-L)  for  L*t 

SlJblP)  * SINP(THETA) 

HARM 

165 

If  ( MT-2 ) 260 . 220  * 220 

HARM 

1«7 

220 

DO  250  l32»  MT 

HARM 

168 

THETA  a T4FTA/?. 

harm 

189 

JSTfeP'2  * JSTEP 

harm 

190 

JSTtP  a JJIf 

harm 

191 

JOlf  ■ JO  I P/2 

S< JiIFjaSJNFlTHETA) 

HARM 

192 

JClaNT-jDJP 

scjld  aC 3sf < theta ) 

HARM 

1*4 

JLAST*NT-JSTEP2 

harm 

1*6 

If  I jLAST* JSTEP) 250.230. 230 

HARM 

1*7 

2 jO 

DO  2«0  J« JSTEP. JLAST. JSTEP 

HAHM 

1*8 

JC»nT-J 

HARM 

1** 

J0»j*JDIF 

HARM 

200 

240 

Si Ji )*S( J)*S< JCi>*S( JUlf i«S(JC) 

HARM 

201 

290 

CONI InUE 

HARM 

202 

260 

If  « IPS>20,1.20 

harm 

203 

EM 

SCOPE 

Ha  00  c3*u 

Ha  00  RJN.200. 60000 


•U.1  OOVERNMCNT  fftMTHM  Of ftC*:*TI  Ml-150/0055 
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